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PREFACE 


This  is  Volume  1  of  a  two  volume  final  report  for  the  contract  "Analysis  of 
Dynamical  Plasma  Interactions  with  High  Voltage  Spacecraft."  The  period  of 
technical  performance  was  31  December  1988  through  31  December  1991.  The 
objectives  of  this  contract  were  to  study  dynamical  plasma  interactions  with 
high  voltage  spacecraft,  to  construct  a  three-dimensional  computer  code, 
DynaPAC,  as  a  workbench  for  such  studies,  and  to  support  the  SPEAR 
program.  Volume  1  is  a  compilation  of  work  done  to  model  high  voltage 
plasma  interactions,  with  application  to  the  SPEAR-II  chamber  tests  and  to 
the  design  of  SPEAR-3.  Volume  2  contains  documentation  for  the  DynaPAC 
code,  as  well  as  for  the  two-dimensional  Gilbert  code. 


1.  INTRODUCTION 


This  volume  contains  excerpts  of  work  done  in  three  general  areas; 

(1)  Dynamic  plasma  problems  done  in  two  dimensions  using  the  Gilbert 
code.  (Chapters  2-3.) 

(2)  Work  performed  in  support  of  the  SPEAR-II  program.  (Chapters  4-7.) 

(3)  Work  performed  in  support  of  the  SPEAR-3  program.  (Chapters  8-10.) 
The  SPEAR-3  work  is  continuing  under  contract  F19628-91-C-0187, 
entitled  "Space  System-Environment  Interactions  Investigation." 

In  addition  to  numerous  technical  meehngs  directly  related  to  this  contract  or 
to  the  SPEAR  program,  presentations  were  made  under  this  contract  at  the 
following  professional  conferences: 

Spacecraft  Charging  Technology  Conference,  Monterey,  CA,  30  October 
through  2  November,  1989.  "Dynamics  of  Spacecraft  Charging  by  Electron 
Beams."  M.  J.  Mandell  and  I.  Katz.  (Conference  proceedings  available.) 

High  Voltage  Workshop,  Port  Hueneme,  CA,  19-21  March,  1991.  "Transient 
Plasma  Current  to  a  High  Voltage  Probe."  M.  J.  Mandell,  T.  T.  Luu,  G.  A. 
Jongeward  and  I.  Katz. 

Physics  Computing  '91, 10-14  June,  San  Jose,  CA.  "DynaPAC  -  A  3-D  Finite 
Element  Plasma  Analysis  Code."  M.  J.  Mandell,  T.  T.  Luu,  and  J.  R.  Lilley 

SOAR  '91,  NASA/JSC,  Houston,  TX,  9-11  July  1991.  "Spacecraft-Plasma 
Interaction  Codes:  NASCAP/GEO,  NASCAP/LEO,  POLAR,  DynaPAC,  and 
EFSAT."  M.  J.  Mandell,  G.  A.  Jongeward,  and  D.  L.  Cooke.  (Conference 
proceedings  available.) 
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Chapters  2  and  3  of  this  report  present  two-dimensional  calculations  of 
generic  plasma  interaction  situations  in  which  dynamic  effects  strongly 
influence  measurable  results.  The  work  in  Chapter  2  was  motivated  by 
observations  by  the  MAIMIK  rocket  of  spacecraft  potentials  far  in  excess  of  the 
energy  of  an  emitted  electron  beam.  This  work  shows  that  the  oscillations  of  a 
beam  emitted  into  an  underdense  plasma  can  cause  this  overcharging. 
Chapter  3  is  an  attempt  to  calculate,  for  fairly  modest  parameters,  the  effects  of 
the  sheath  instability  caused  by  orbital  motion  of  a  positive  body.  It  shows 
that  substantial  potential  oscillations  occur,  and  that  the  electron  current  to 
the  probe  is  enhanced. 

Chapters  4-7  present  work  performed  in  support  of  SPEAR-Il.  Chapter  4  deals 
with  the  generic  problem  of  plasma  response  to  a  sudden  application  of 
negative  voltage  to  a  probe.  The  calculation  was  done  first  with  Gilbert,  then 
duplicated  with  DynaPAC.  It  showed  that  it  takes  a  long  time  to  establish  an 
equilibrium  sheath,  and  that  during  this  time  ion  currents  to  the  probe  can  be 
enhanced  by  as  much  as  an  order  of  magnitude.  Chapter  5  shows  the 
application  of  these  ideas  to  the  sheath  of  the  SPEAR-U  high  voltage  probe, 
explaining  the  observed  difference  in  probe  reading  between  vacuum  and 
plasma  conditions.  Actual  DynaPAC  calculations  (performed  post-test  during 
the  fall  of  1990)  of  the  actual  plasma  currents  to  the  SPEAR-II  high  voltage 
components  appear  in  Chapter  6.  (It  is  worth  noting  that  these  calculations 
proceeded  at  a  rate  of  one  one-microsecond  timestep  per  day  on  a  Sun 
Microsystems  SPARCStation  I.  Subsequently,  DynaPAC  recoding  led  to  at 
least  a  factor  of  three  speed  increase,  and  transfer  to  an  SGI  Iris  Indigo  to 
another  factor  of  three.  Were  we  to  repeat  the  calculation,  we  should  be  able 
to  do  several  timesteps  per  day.)  The  SPEAR-II  calculations  are  tied  together 
and  summarized  in  S-Cubed's  contribution  (partially  supported  by  other 
contracts)  to  the  SPEAR-II  Products  Report,  which  is  included  as  Chapter  7. 
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Finally,  Chapters  8-10  show  some  early  work  on  the  SPEAR-3  modeling. 
Chapter  8  describes  a  one-dimensional  model,  since  incorporated  into  the 
EPSAT  code,  which  provided  early  estimates  for  the  neutral  density  required 
for  breakdown  of  a  negative  sheath.  Chapters  9  and  10  describe  early  work  on 
the  rocket  floating  potentials,  collected  currents,  and  current  distribution  to 
the  particle  detectors.  With  the  benefit  of  this  work,  we  were  able  to  improve 
our  computational  models  and  update  our  geometrical  models  to  obtain 
excellent  results  under  the  new  contract. 
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2.  DYK  AMiCS  OF  SPACECRAFl  CHARGING  BY  ELECTRON  BEAMS 


This  work  was  presented  at  the  Spacecraft  Charging  Technology  Conference, 
Monterey,  CA,  30  October  through  2  November,  1989.  It  was  published  in  the 
conference  proceedings.  It  also  appeared  in  the  Interim  Report  (30  September 
1990)  for  this  contract. 
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Dynamics  of  Spacecraft  Charging  by  Electron  Beams 


M,  J.  Mandell  and  I.  Katz 

S-CUBED  Division  of  Maxwell  Laboratories,  La  Jolla.  California 

When  a  spacecraft  or  rocket  emits  an  electron  beam  into  an  undcrdcnsc  plasma,  it  can 
charge  to  potentials  in  excess  of  the  beam  energy.  We  show  calculations  in  which  8  kcV 
beams  are  emitted  along  and  across  the  earth's  magnetic  field,  with  parameters  appropnate 
to  the  MAn^OK  rockeL  As  was  observed  on  M.AIMIK,  the  spacecraft  charged  to 
potentials  in  excess  of  the  beam  potential  due  to  energization  of  beam  electrons  by  beam- 
generated  electrostatic  oscillations.  This  is  in  contrast  to  the  low  levels  of  charging  often 
seen  in  denser  environments,  where  higher  plasma  currents,  coupled  with  ionization  of 
neutrals,  hold  spacecraft  potentials  below  a  few  hundred  volts. 

The  beam  structure  and  the  spacecraft  potential  oscillate  at  a  frequency  corresponding 
to  the  mean  lifetime  of  the  beam  electrons.  These  oscillations  energize  the  beam 
electrons,  and  also  pump  energy  into  the  ambient  plasma,  which  exhibits  lower  frequency 
oscillations.  The  peak  spacecraft  potential  is  over  1  kV  in  excess  of  the  beam  energy. 

For  the  cross-field  case,  the  oscillation  frequency  is  proportional  to  the  beam  current  for 
sufficiently  intense  beams. 

Following  beam  tum-off,  there  is  an  immediate  return  of  unscattercd  beam  electrons 
and  a  longer  term  dissipation  of  scattered  beam  electrons.  Analytic  estimates  are  presented 
for  the  decay  and  overshoot  of  the  spacecraft  potential. 


INTRODUCTION 

A  number  of  electron  beam  experiments  have  measured  results  that  apparently  violate 
conservation  of  energy.  The  SEPAC  experiment  (Reasoner  et  al.,  1984)  measured  a 
spectrum  of  returning  electrons  extending  to  energies  well  above  that  of  the  emitted  beam, 
and  the  MAIMIK  rocket  (Maehlum  et  al.,  1988;  Denig,  Maehlum  and  Svencs,  this 
conference)  was  charged  by  an  8  keV  electron  beam  to  potentials  as  high  as  14  keV.  Katz 
et  al.  (1986)  performed  a  planar  calculation  that  showed  that  oscillations  of  the  beam 
electrons  led  to  the  spectral  broadening  seen  in  the  SEPAC  experiment.  In  this  paper,  we 
show  2-dimensional  calculations,  with  parameters  appropriate  to  MAIMIK,  illustrating  that 
space  charge  oscillations  associated  with  the  electron  beam  lead  to  rocket  potentials  in 
excess  of  the  beam  energy. 

The  electron  beam  on  MAIMIK  was  directed  nearly  across  the  earth’s  magnetic  field. 
The  geometry  of  an  intense  beam  directed  across  a  magnetic  field  is  shown  in  Figure  la.  If 
the  spacecraft  is  near  the  beam  energy,  the  beam  electrons  will  be  slowest  at  their  farthest 
excursion  from  the  spacecraft  and  will  form  a  space  charge  maximum.  This  space  charge 
maximum  breaks  the  azimuthal  symmetry  of  the  problem,  so  that  electrons  may  be  scattered 
from  their  original  gyro-orbits  and  leave  the  vicinity  of  the  spacecrafL  Unfortunately,  this 
is  a  truly  3-dimensional  situation  and,  therefore,  very  difficult  to  model. 

Two  beam  configurations  that  can  be  modeled  in  2-dimensional  axisymmetric  geometry 
are  shown  in  Figures  lb  and  Ic.  Figure  lb  shows  a  beam  directed  along  the  magnetic 
field.  An  intense  virtual  cathode  is  formed.  Most  of  the  beam  electrons  move  outward 
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frcni  the  virtual  cathcxle  and  are  attracted  back  to  the  spacecraft,  which  they  impact  from 
behind.  A  few  of  the  beam  electrons  continue  along  the  magnetic  field  line  and  escape  the 
vicinity  of  the  spacecraft.  Figure  Ic  illustrates  an  "equatorial''  beam,  in  which  electrons  are 
directed  across  the  magnetic  field  from  tne  entire  spacecraft  circumference.  A  ring  of 
maximum  space  charge  appears  around  the  spacecraft.  Because  the  space  charge  maximum 
is  a  full  ring  rather  than  a  localized  region,  we  expect  it  to  be  less  intense  than  would  be  the 
case  for  a  physical  beam  of  the  same  current.  Also,  the  space  charge  ma.ximum  is  less 
effective  at  scattering  because  it  does  not  break  azimuthal  symmetry. 

The  calculations  described  here  were  performed  using  an  S-CUBED-deve!oped,  finite- 
element,  electrostatic  particie-in-cell  code  named  Gilbert.  Gilbert  is  a  flexible,  multi¬ 
purpose  code  with  many  special  features.  For  these  calculations,  the  space  around  the 
spacecraft  was  gridded  with  biquadratic  elements  of  variable  resolution  to  a  distance  of  ten 
meters.  Each  element  represents  a  volume  of  space  corresponding  to  its  area  revolved  in  a 
full  circle  about  the  symmetry  axis.  For  the  "equatorial"  beam,  the  calculation  took 
advantage  of  mirror  symmetry  about  the  equatorial  plane.  The  computational  grids  for  the 
two  cases  are  shown  in  Figure  2. 

PHYSICAL  PARAMETERS 

Table  1  shows  the  physical  parameters  used  in  the  calculation.  The  parameters  listed  in 
the  "current  collection"  category  are  calculated  for  a  spherical  spacecraft  at  the  beam 
potential  of  8,000  volts. 


Table  1.  Ffroblem  Parameters 


Geometrical  Parameters 
Inner  (spacecraft)  Radius  0.3  m 

Outer  Radius  10.0  m 

Magnetic  Field  0.4  gauss 

Beam  Parameters 

Beam  Current  0.16  amperes 

Beam  Energy  8000  eV 


Plasma  Parameters 
Electron/Ion  Density 
Electron  Temperature 
Ion  Mass 

COpe 

“ce 

J.h 


3  X  10^  m-3 
leV 
16  amu 

13.6cm 
3.1  X  10^  sec'^ 
7.0  X  10®  sec' 
8.0  X  10  ®  A-m  '^ 


Current  Collection 
Langmuir-Blodgett  Radius 
Parker-Murphy  Radius 
Parker-Murphy  Current 
Probe  Charge 


10  m 
1.54  m 

0.0025  amperes 
2.7  X  10  coulombs 
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The  plasma  surrounding  the  spacecraft  is  "underdense  ',  in  the  sense  that  the  Parker 
Murphy  bound  (Parker  and  Murphy,  1967)  on  the  plasma  return  current  is  only  a  few 
percent  of  the  beam  current.  Also,  the  electron  plasma  frequency  is  below  the  electron 
cyclotron  frequency  and,  as  we  shall  see,  well  below  the  oscillation  frequency  of  the  beam. 

The  radius  of  the  computational  space  is  taken  equal  to  the  radius  of  a  space-charge- 
limited  spherical  sheath  in  this  plasma,  which  is  far  larger  than  the  Parker-Murphy  radius 
for  current  collections.  The  effect  of  the  plasma  external  to  this  radius  is  represented  by  a 
zero  potential  condidon  on  this  boundary.  In  retrospect,  this  approximation  appeon 
adequate  to  represent  the  beam  dynamics,  the  ion  dynamics,  and  the  collecdon  of  ambient 
electrons  but  omits  long-range  transient  effects  on  the  ambient  electron  density. 

CALCULATIO.N  FOR  BEAM  ALONG  FIELD 


The  calculation  begins  with  the  grid  of  Figure  2a  filled  with  electron  and  ion 
macropardcles  (Tigure  3a).  The  beam  iS  projened  along  the  magnetic  field  (Figure  3b)  and 
initially  exit.s  the  computational  space.  The  negative  beam  charge  and  the  positive  charge 
left  behind  on  the  spacecraft  produce  a  dipole  potential  (Figure  3c)  that  expels  ambient 
electrons  from  the  beam  region  (Figure  3d} 

The  spacecraft  reaches  beam  potential  about  2.5  microseconds  after  beam  lum-on 
(Figure  4a),  and  the  potential  exhibits  persistent  oscillations  (at  2  x  10*  Hz)  for  the  duration 
of  the  calculation  (Figure  4b).  During  this  time,  the  spacecraft  remains  always  above  the 
beam  potential  and  has  a  peak  potential  of  about  9,100  volts. 

The  beam  conformation  oscillates  along  with  the  spacecraft  potential.  Figure  5a  shows 
the  beam  conformation  when  the  spacecraft  potential  is  fairly  high.  The  bulk  of  the  beam 
electrons  are  far  from  the  spacecraft,  having  been  emitted  when  the  potential  was  low. 
Figure  5b  shows  beam  conformation  at  a  fairly  low  potential.  In  this  case,  the  bulk  of  the 
beam  electrons  are  close  to  the  spacecraft,  having  been  emitted  at  high  potential.  However, 
a  pulse  of  energized  electrons  can  be  seen  escaping  along  the  field  line. 

Figures  6a  and  6b  show  the  electrostatic  potential  structure  about  the  spacecraft.  The 
sheath  is  elongated  along  the  magnetic  field  due  to  quasi-trapping  of  ambient  electrons  that 
cannot  be  collected.  The  contours  are  distorted  along  the  axis  by  the  beam  electron  space 
charge  (Figure  6a).  At  times,  a  negative  potential  well  forms  in  the  cross-field  region 
(Figure  6b). 

Figure  7a  shows  the  amount  of  beam  electron  charge  in  the  computational  space,  and 
Figure  7b  shows  the  time  dependence  of  the  dipole  moment  of  the  beam  electrons,  defined 
as 


Dipole  Moment  =  I  p  z  d-^r 


I 


where  p  is  the  charge  density  of  beam  electrons.  Since  the  dipole  moment  is  oscillating  at 
2  MHz,  we  expect  to  see  strong  electromagnetic  radiation  at  this  frequency. 

Figures  8abc  show  spectral  analysis  of  some  of  the  oscillating  quantities.  1  he  potential 
(Figure  8a)  shows  a  sharp  peak  at  2  MHz,  with  well  -defined  second  and  third  harmonics. 
The  beam  dipole  moment  (Figure  8b)  shows  a  sharp  peak  at  2  \THz,  as  well  as  a  broad 
peak  at  the  ambient  electron  plasma  frequency  (0.5  MHz).  The  dipole  moment  of  the 
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amb:cni  :  h!.:’!rtr  •'^c’s  sf.o'.vs  iit:;'!; 


h'cnvT  of  {he  1  N!H/  osc-ilat-ons  but  h;s'v  a 


strop;’,  Hro:\i1  peak  near  the  antb;ent  plasriu  frequency 

The  beam  charite  (Figure  7a t  vanes  from  t).()4  to  0  08  jaeou!  iJividing  this  by  the 
beam  current  ot'O  16  amperes  shows  that  beam  parocic  iifcumes  fall  mainly  in  the  range 
0.25  0.5  pscc.  Noting  also  that  the  upstrokes  m  Figure  1 1  strong  return  current)  are 
steeper  than  the  downstiokes  (sieailily  emitted  current;  leads  to  tiic  following  intcrpretanon 
Most  of  the  beam  electrons  returii  to  the  spacecraft  m  a  shon  burst,  F'cctr(.ins  emitted 
during  the  return  burst  see  relatively  weak  retarding  fields,  travel  far  from  the  spacecraft, 
and  have  lifeumes  of  approximately  a  full  oscillation  penod  Flccnrons  emitted  while  the 
potential  is  nsing  see  stronger  retarding  fields,  travel  smaller  excu'‘s!ons.  and  return  lo  the 
spacecraft  at  the  same  ume  as  the  long-lived  electrons,  producing  the  return  cunent 
bunching.  This  bunching  is  similar  to  that  seen  by  Katz  ei  al  ( 1086) 

Figure  9  shows  the  ion  positions  at  the  conclusion  of  the  calculation  (20  pseci.  The 
ions  have  clearetd  a  region  of  about  two  meters  .iround  die  spacecraft.  Thus  the  run  time  of 
this  calculation  is  tcxi  shon  to  approach  the  fonnation  of  an  equilibrium  sheath  ‘ 


(  ai.c:ti..\tion  for  bk.\m  ac  ross  fift.d 


A  similar  calculation  w  a.s  pen'omied  for  an  equatonal  beam  directed  across  the  magneuc 
field  in  the  grid  of  Figure  2b.  (Note  that  in  this  calculation  we  have  taken  advantage  of 
mirror  symmetry'  aboc’  the  z  -  0  plane.  \'aiues  for  current  and  chtirge  will  be  quoted  at 
double  the  computed  values,  so  that  they  are  charactensne  of  a  complete  sphere.) 

ITie  beam  electrons  at  the  conclusion  of  the  calculation  are  shown  in  Figure  10.  The 
figure  shows  a  main  stream  of  electrons  emitted  from  the  spacecraft,  slowed  by  the  electric 
field  and  ntmed  by  the  magnetic  field,  then  returning  to  the  spacecraft.  In  addioon,  there  is 
a  column  of  beam  electrons  extena.ng  in  the  magnetic  field  direction  The  mechanism  for 
populating  this  column  is  that  a  beam  electron  in  iLs  midal  orbit  gains  enough  momentum 
along  the  magnenc  field  to  miss  the  spacecraft  on  its  first  return  passage,  and  while  passing 
near  the  spacecraft  receives  a  substantial  impulse  along  the  field.  Electrons  leave  the 
column  by  either  impacung  the  spacecraft  or  escaping  the  gnd. 

Figures  1  lab  show  the  time  history  of  the  spacecraft  potential.  Four  different  current 
values  were  used.  From  Figure  I  lb  we  see  that,  while  doubling  the  current  led  to  some 
increase  in  the  mean  potential  and  its  oscilianon  amplitude,  the  main  effect  is  to  double  the 
frequency. 

Figures  12ab  show  the  potential  contours  about  the  spacecraft  at  two  different  dmes. 

As  in  the  previous  case,  the  sheath  is  elongated  along  the  field  line,  and  a  negative  potential 
*  i.>.  ' .  lAimc'iinv  in  the  cross-field  region. 

Figure  13  shows  u.e  ambient  electron  macroparecles.  The  ambient  electron 
population  was  maintiuned  by  generating  the  plasma  thennal  current  at  the  spheni.al 
problem  boundary  out  to  a  radius  of  eight  meters.  A  low  density  region  is  seen  to  extend 
along  the  field  line  from  the  spacecraft;  electrons  in  this  region  are  allowed  by  the  theory 
of  Parker  and  Murphy  ( 1967)  to  be  collected  by  the  spacecraft.  A  high-density  cloud  of 
electrons  is  seen  in  the  cToss-field  region,  as  these  electrons  cannot  be  collected  by  the 
spacecraft  and  have  low  probability  of  escaping  the  gnd.  This  ambient  electron  charge 
density  structure  is  the  cause  of  the  elongated  potential  contours  shown  in  Figures  12ab 
and  6ab. 
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dv/di  =  a(/c)  (v/v^i* 


where  is  the  Parker-Murphy  current  at  poientia]  Vg  (-  0,0025  amperes  at  8,000  volts) 
and  C  is  the  spacecraft  capacitance  (3.4  x  10  farati).  The  solution  is 

tj-t,  =-2.4x  lO  ^’  V,^^) 

which  gives  a  discharge  time  of  about  200  psec.  Since  the  thermal  ion  speed  is  only  a  few 
thousand  m-sec  ^  it  will  take  about  a  millisecond  for  the  ions  to  repopulatc  a  sheath  a  few 
meters  in  radius.  Assuming  a  thermal  electron  distribution,  the  spacecraft  will  charge 
negatively  during  this  time  according  to 

dV/dt  =  47ta2  4  e'^^/C 


whose  solution  is 

V/0  =  -  /n(l  +  t/x) 
t  =  C0/(4tta2|4  I) 

For  our  parameters,  x  is  about  0.4  psec,  and  we  expect  an  overshoot  to  about  -8  volts. 
However,  if  the  electron  distribution  has  an  elevated  thermal  tail,  as  is  likely  due  to 
turbulence  associated  with  the  nonuniform,  nonequilibrium  ion  distribution,  the  negative 
overshoot  will  be  much  greater. 


CONCLUSIONS 

Electron  beams  emitted  from  spacecraft  in  the  ionosphere  exhibit  complex  behavior. 

We  have  analyzed  here  the  case  of  a  beam  emitted  into  an  underdense  plasma,  with 
parameters  appropriate  to  the  MAIMIK  rocket.  The  beam  was  emitted  both  along  the 
magnetic  field  and  across  the  field  in  an  "equatorial"  fashion. 

The  beam-emitting  system  exhibits  oscillations  at  a  few  megahertz.  These  oscillations 
are  associated  with  bunched  return  of  the  beam  electrons  and  cause  electron  energization  so 
that  the  spacecraft  can  achieve  potentials  in  excess  of  the  beam  energy.  At  least  for  the 
field-aligned  beam  case,  these  oscillations  are  dipolar  in  character  and  should  be  observable 
as  electromagnetic  radiation.  For  sufficiently  intense  beams,  the  oscillation  frequency  is 
proportional  to  the  beam  current. 

The  dynamics  of  the  ambient  plasma  is  largely  independent  of  the  beam  dynamics.  The 
ambient  electrons  show  a  broad  peak  at  their  own  plasma  frequency  and  form  an  elongated 
sheath. 

Even  in  this  very  underdense  plasma,  the  relaxation  of  spacecraft  potential  following 
beam  turn-off  is  rapid  compared  with  the  time  for  ions  to  thermally  fill  in  the  sheath.  If  the 
nonuniform  plasma  that  exists  during  this  time  causes  electron  heating,  a  substantial 
negative  overshoot  of  the  spacecraft  potential  can  occur. 

Acknowledgement.  This  woric  was  supported  by  the  Geophysics  Laboratory,  Hanscom  Air  Force  Base, 
Massachusetts,  under  contract  F19628-89-C-0032. 
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Tabie  2  jp.ow  s  ihe  range  of  beani  charge,  range  of  beam  cieccron  lifeariie,  and 
oscillation  penod  for  the  dtffere.nt  values  of  errassion  current.  As  with  the  ficld-aiigned 
beam,  the  oscillauon  period  is  near  the  maxunum  beam  particle  lifedme.  Except  for  the 
lowest  current  value,  the  beam  charge  increases  only  slightly  when  die  current  is  doubled. 
This  is  because  the  maximum  beam  charge  is  approaching  the  charge  on  the  sphere.  When 
the  current  is  doubled,  a  modest  increase  in  mean  spacecraft  potential  is  sufficient  to  reduce 
the  particle  excursion  distance  and  cut  the  beam  panicle  lifenme  by  half.  This  point  is 
further  illustrated  by  Figure  14.  which  plots  the  quantity 

I  prd’r 

The  average  and  oscillation  amplitude  of  this  quandty  varies  slow  ly  as  the  sheath  is  being 
formed,  but  shows  no  abrupt  changes  as  the  current  is  altered,  leading  to  the  conclusion 
that  the  system  rapidly  adjusts  so  that  the  beam  charge  umes  its  excursion  distance  is 
independent  of  current. 

Table  2.  Beam  Charge.  Lifetime,  and  Oscilladon  Period  for  Equatorial  Beam 


Current 

Charge 

lifetime 

Period 

f.^rriDeresl 

(ucoul) 

(use.c) 

0.16 

0.07  -  0.1 

0.44  -  0.62 

0.62 

0.32 

0.12  -  0.2 

0.38  -  0.62 

0.47 

0.64 

0.13  -  0.21 

0.20  -  0  33 

0.28 

1.28 

0.14  -  0.25 

0.11  -  0.19 

0.18 

BEAM  TURN-OFF 

Several  papers  at  this  conference  discussed  beam  turn-off,  showing  that  the  spacecraft 
poiendal  tends  to  ovenhoot  and  achieve  a  negative  value.  This  led  us  to  investigate  the 
behavior  of  the  equatorial  beam  system  after  turn-off  of  the  0. 160  ampere  beam. 

Figure  1 5  shows  the  behavior  of  the  spacecraft  potential  and  the  beam  charge  following 
beam  turn-off.  Most  of  the  beam  charge  is  promptly  collected,  dropping  the  spacecraft 
poiendal  to  7,200  volts.  Approximately  0.03  pcoul  of  charge  remains  in  the  field-aligned 
column  of  scattered  beam  electrons,  which  decays  with  a  time  constant  of  2.1  p5cc.  About 
three-quarters  of  this  depopulation  rate  results  from  electrons  escaping  the  grid,  and  about 
one-quarter  from  electrons  recaptured  by  the  spacecraft.  (When  the  spacecraft  is  at  elevated 
potential,  the  escape  rate  must  balance  the  Parker-Murphy  collected  current,  leading  to  a 
12  p.sec  time  constant  for  escape.) 

If  the  spacecraft  discharges  by  collecting  the  Parker-Muiphy  current,  its  potential  will 
follow  the  equation 
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Fig,  1 ,  Three  beam  configiiracions  from  a  spherical  spacecraft:  (a)  physical  tx-am  across  magnetic  fielil,  (h) 
physical  beam  along  magnetic  field;  (c)  "equatorial"  beam  across  magnetic  field 
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Fig.  2.  Axisymmetic  gnds  of  bu|uadratic  elemcnis  used  in  these  calculations:  (a)  gnd  used  for  calculation  of 
beam  along  magnetic  field  (fig.  la);  (b)  gnd  used  for  calculation  of  "equatorial '  beam  (fig.  lb)  using  mirror  plane. 
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Figure  3.  (a)  Initial  electron  or  ion  inacroparticles  in  gnd  of  fig.  2a;  (b)  beam  electron  macroparticles 
p.sec;  (c)  dipolar  potential  at  0.04  psec,  with  contours  marked  A,  B,  and  C  at  -50,  -2,  and  +5  voll.s. 
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Fig.  4.  Time  dependence  of  spacecraft  potential  for  field-aligned  beam:  (a)  initial  transient;  (b)  persistent 
oscillations. 
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Fig.  5.  Beam  conformation  (a)  at  high  poienlial 
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Fig.  6.  Electrosiatic  potential  structure  (logarithmically  spaccti  co 
(b)  with  negative  potential  well  in  the  cross-field  region. 
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Fig.  7.  Time  dependence  of  (a)  beam  electron  charge;  and  (b)  beam  electron  dipole  moment. 
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Fig.  12.  Electrostatic  potential  structure  (loganthmically  spaced  contours) 
(b)  with  negative  potential  well  in  the  cross-field  region 
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Fig.  13.  Ambient  e! 


quantity  I  p  r  d  r  for  the  beam  electrons,  showing  the  lack  of  dependence  on  stn  ng  current 


Fig.  15.  Behavior  following  beam  turn-off:  (a)  spacecraft  potential,  (b)  beam  electron  charge. 


3.  POSITIVE  ORBITAL  PROBE  CALCULATION 


This  work  appeared  in  the  Quarterly  Report  for  1  January  through  31  March 
1991. 
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3.  POSITIVE  ORBITAL  PROBE  CALCULATION 


3.1  RATIONALE 

We  are  accustomed  to  perlormtng  “sheath"  calculations  of  currents  to  spacecraft. 
The  “sheath"  is  defined  as  the  region  (or  boundary  of  the  region)  from  which  the 
repelled  species  is  excluded.  A  common  and  apparently  benign  case  for  which 
this  entire  concept  makes  no  sense  is  an  orbiting  ionospheric  probe  biased  at, 
say.  +10  volts.  The  repelled  species  (O-^  ions)  will  penetrate  the  potential 
structure  to  the  5-volt  level.  In  order  to  understand  the  physics  of  this  situation, 
we  attempted  to  use  Gilbert  to  perform  a  full-PIC  simulation  of  this  situation. 


3.2  CALCULATION  DESCRIPTION 


Calculation  parameters  were  chosen  so  that  (1)  the  problem  is  realistic;  (2)  the 
zone  size  would  not  exceed  a  few  Debye  lengths;  (3)  the  range  of  electron 
energies  would  be  fairly  tolerable;  and  (4)  time-stepping  on  the  electron 
timescale  would  eventually  lead  to  correct  ion  motion.  The  calculation 
parameters  (along  with  several  derived  quantities)  are  shown  in  Table  3.1. 


Table  3.1  Parameters  for  the  positive  orbital  probe  calculation. 


Probe  Radius 

0.1  m 

Probe  Potential 

10  V 

S/C  Velocity 

7500  m/sec 

Ram  Energy  (0+) 

5  eV 

Magnetic  Field 

None 

Plasma  Density 

10'*''  m'3 

cope 

1.78  X  107 

log  10  Wpe 

7.25 

Electron  Temperature 

0.1  eV 

1  eV 

Debye  Length 

7.44  X  10'3  m 

2.35  X  10-2  m 

Elec.  Th.  Curr. 

8.46  X  lO  "^  A/m^ 

2.68  X  10-3  Aym2 

L.-B.  Radius 

0.26  m 

0.22  m 

L.-B.  Current 

0,72  mA 

1.63  mA 

27 


r 


VVe  have  included  derived  parameters  for  1  eV  as  well  as  0. 1  eV  electron 
temperature  because  it  can  be  argued  that  the  boundary  conditions  had  the 
effect  of  raising  the  effective  electron  temperature.  However,  the  calculated 
electron  current  to  the  probe  was  3.84  mA,  which  is  more  than  double  the  probe 
current  in  a  stationary  1  eV  plasma. 

Figure  3.1  shows  the  axisymmetric  computational  grid  for  the  problem. 
Conductor  1  represents  the  probe.  Conductor  2  (0.5  m  upstream  of  probe 
center)  emits  electrons  and  ram  ions.  Conductor  3  (0.5  m  radially  from  probe 
center)  and  conductor  4  (0.3  m  downstream  from  probe  center)  emit  electrons 
only.  The  potential  at  conductors  2,  3.  and  4  are  held  fixed  at  0  volts,  while 
conductor  1  is  held  fixed  at  +10  volts. 


Figure  3.1  Axisymmetric  Computational  grid  for  positive  potential  orbital 
probe.  Grid  is  50  cm  in  radius  and  80  cm  in  height. 
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Even  though  boundaries  2  and  3  are  well  beyond  the  Langmuir-Blodgett  sheath 
distance,  it  was  necessary  to  emit  more  than  the  plasma  thermal  current  to  avoid 
a  sharp  potential  rise  at  the  boundary.  The  emission  used  followed  the  relation 

J  =  2.33  X  10-6  v3/2/d2  +  2.12  X  10^  ne 


where  n  is  the  ambient  ion  density  and  V  is  the  potential  at  the  zone  boundary 
at  distance  d.  This  equation  comes  from  treating  a  constant  ion  density  as  a 
perturbation  to  the  Child-Langmuir  problem. 

The  calculation  proceeded  by  first  running  long  enough  for  the  ions  to  fill  the 
grid  and  establish  an  apparent  steady  state,  then  resetting  the  clock  and 
continuing  to  run  and  take  data. 


3.3  CALCULATION  RESULTS 

Figure  3.2  shows  a  snapshot  of  the  ion  positions  toward  the  end  of  the  calculation. 
The  stagnation  line  is  well-defined  from  about  30  to  90  degrees.  Below  30 
degrees,  a  high  amplitude  of  electrostatic  noise  blurs  the  stagnation  line.  Because 
(as  we  shall  see)  the  stagnation  line  is  a  potential  maximum,  it  is  blurred  beyond 
90  degrees  because  ion  trajectories  cannot  stably  follow  the  potential  maximum. 


Figure  3.2  Ion  macroparticle  positions  for  positive  orbital  probe. 
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Figures  3.3  (a  and  b)  show  potential  contours  at  two  slightly  different  times.  The 
potential  maximum  along  the  stagnation  line  (which  shows  up  clearly  m 
Gouraud-shaded  potential  plots)  is  seen  as  an  extension  of  the  +1  volt  contour 
toward  an  island  of  +1  volt  potential  at  the  upper  right.  Part  of  the  stagnation  line 
is  held  below  +1  volt  due  to  the  influence  of  the  electron-rich  wake  region. 
Figures  3.4  (a  and  b)  plot  potentials  along  radial  lines  at  z  =  0.0  and  z  =  0  2 
(relative  to  the  probe  center),  showing  more  clearly  the  potential  maximum  at 
the  stagnation  line,  occurring  at  r  =  0.21  and  r  =  0  31  respectively. 


Marksd  Lavals 


A  -1.00 
B  3.00 
C  7  00 
D  11  00 


A  2  00 
B  200 
C  6  00 
D  10  00 


(a) 


(b) 


Figure  3.3  Potential  contours  around  positive  orbital  probe  for  two  slightly 
different  times.  Spacecraft  velocity  is  downward,  or.  alternatively, 
plasma  flow  is  upward. 
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Potential  Potential 


Both  figures  3.3  (a  and  b)  show  regions  of  negative  potential  m  the  upstream 
direction.  These  regions  are  formed  on  the  axis  and  grow  outward.  The 
timescale  for  this  process  is  the  electron  plasma  frequency.  Figures  3.5  (a,  b.  c, 
and  d)  show  the  time  and  frequency  domain  behavior  of  the  potential  at  node 
27  (3.3  cm  upwind  of  the  sphere  surface)  and  at  node  1 22  ( 1 5  cm  upwind  of  the 
sphere  surface).  At  node  27.  the  noise  level  is  tremendous,  with  a  clear  peak  at 
the  plasma  frequency.  Further  from  the  sphere,  at  node  122,  the  amplitude  of 
the  oscillations  is  much  lower,  and  the  frequency  behavior  much  cleaner. 


12. 


-2. 


-4.4 - . - - - , - ^ - - - - - 1 
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Time 


(a) 

Figure  3.5  Potential  fluctuations  in  time  and  frequency  domain;  (a)  3.3  cm 
upstream,  time  domain;  (b)  3.3  cm  upstream,  frequency  domain; 
(c)  15  cm  upstream,  time  domain;  (d)  15  cm  upstream,  frequency 
domain. 
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Potential  at  Node  1 22 


34 


!I* 


3.4  DOUBLE-LAYER  MODEL 


Noting  that  we  have  a  h.gh  current  to  the  sphere,  we  are  motivateO  to  tne 
folfowing  calculation.  Suppose  we  nave  (tor  smtpliaiy)  a  planar  sheath  with 
ions  penetrating  to  the  ievei  t-  Regarding  the  plasma  as  Deing  a  coio  electron 
source  at  the  O-potential  level,  how  much  electron  current  is  requireo  to  make  to 
zero-to-five  (i.e.,  r)  volt  region  a  stable  double  layer,  and  how  thick  is  the  layert> 

The  double  layer  has  incident  and  reflected  lon  currents  and  transmitted 
electron  currents,  so  that  Poisson's  equation  for  the  layer  region  is 

-t'o  0”  =  po  {  2  lL-(f  -  0)]  ■  2  -  (B  4;i0)f  ^  } 

where  po  is  the  ambient  plasma  density  times  the  unit  charge,  and  U  ito  be 
determined)  is  the  temperature  required  of  the  plasma  for  ns  one-sided  thermal 
current  to  equal  the  m.tnimum  current  for  a  stable  double  layer.  Multiplying  botn 
sides  by  o'  enables  us  to  integrate  this  equation  once.  If  we  then  apply  the 
boundary  condition  that  o'  =  0  at  o  =  0  and  at  o  =  t'  we  find  that  0  =  I6m-  and 

(0')2  =  (2po-'‘‘-'o)  {  4(e{h  -  0)1’  2  +  (eo/;t)h2 . . 

For  a  ram  energy  of  5  eV,  a  stable  double  layer  at  the  sheath  requires  the 
electron  thermal  current  of  a  250  eV  plasma'!!  Since  the  actual  temperature  is 
0.1  eV,  we  can  only  achieve  such  a  current  with  a  "convergence  factor"  of  50, 
This  explains  why  an  extended  sheath  (or  presheath)  must  exist  for  ram 
electron  collection.  For  the  plasma  density  used  above,  we  can  also  show  that 
the  thickness  of  a  stable  double  layer  would  have  a  tower  bound  of  about  3  cm 
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4.  TRANSIENTT  SHEATH  DYNAMICS 


This  work  appeared  in  the  Interim  Report  (30  September  1990)  for  this 
contract.  It  had  previously  been  presented  at  a  SPEAR-II  review,  and  was 
subsequently  presented  at  the  High  Voltage  Workshop,  Port  Hueneme,  CA, 
19-21  March,  1991. 
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4.  TRANSIENT  SHEATH  DYNAMICS 


We  have  long  suspected  that  SPEAR  II  would  not  establish  an  equilibrium  sheath 
because  of  the  tong  time  required  for  ions  to  clear  the  sheath.  The  equilibrium 
sheath  radius  for  SPEAR  II  is  about  two  meters,  while  the  most  optimistic  ion 
speed  estimate  (100  keV  ions)  is  1  meter/psec.  This  observation  motivated  us  to 
study  transient  sheath  dynamics  using  the  2-dimensional  Gilbert  code.  These 
results,  with  a  videotape  showing  particle  and  sheath  dynamics,  were  presented 
at  the  SPEAR  11  review  held  at  Chandler,  Arizona,  March  14,  1990. 


4.1  PHYSICAL  AND  COMPUTATIONAL  PARAMETERS 


The  parameters  for  this  calculation  are  shown  in  Table  4.1. 


Table  4.1  Parameters  for  transient  sheath  calculation. 


Sphere  radius 
Boundary  radius 
Sphere  potential 
Plasma  density 
Plasma  temperature 
Debye  length 
Ion  species 
ion  thermal  current 
Timestep 


0.3  m 
10m 

-10®  (1  -  e-‘'^  5  <  io-«) 
10”  m'® 

0  1  eV 
.00743  m 
0+ 

4.94  X  10'®  A-m'2 
2.5  X  10'®  sec 


The  calculation  was  done  for  a  0.3-meter  sphere  in  a  10-meter  boundary. 
Plasma  conditions  were  appropriate  to  the  SPEAR  II  flight.  The  voltage  risetime 
was  chcsen  to  be  appropriate  to  the  3  and  10  psec  SPEAR  II  pulses.  The 
equilibrium  (Langmuir-Blodgett)  solution  to  this  problem  is  a  sheath  radius  of 
15  meters  and  a  current  of  14  milliamperes. 

The  calculation  was  started  with  ion  macroparticles  representative  of  a  uniform 
charge  density.  Electrons  were  taken  to  be  barometric  for  negative  potentials 
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and  to  increase  linearly  for  positive  potentials.  No  "charge  stabilization"  was 
used.  Because  of  the  nonlinear  charge  density,  several  potential  iterations 
were  needed  for  each  0.025  psec  timestep  during  the  early  parts  of  the 
calculation,  during  which  the  applied  potential  was  rising  rapidly.  The  ion 
charge  absorbed  by  the  sphere  was  accumulated,  allowing  the  current  as  a 
function  of  time  to  be  deduced.  The  sheath  radius  was  determined  by  periodic 
visual  inspection  of  the  potential  plots. 


4.2  RESULTS  AND  DISCUSSION 

The  results  of  the  calculation  (current  and  sheath  radius  vs.  time)  are  shown  in 
table  4.2.  The  sheath  expands  rapidly  to  four  meters  as  the  potential  is  applied, 
then  grows  slowly  due  to  erosion.  The  current  peaks  at  about  1 5  psec  at  a 
value  to  times  the  equilibrium  current,  and  then  drops  slowly  thereafter. 

Table  4.2  Results  of  the  transient  sheath  calculation. 


Time 

Potential 

Current 

Sheath  Radius 

[psec] 

[kV] 

[mA] 

[m] 

0.5 

-28 

4 

1.5 

1.0 

-49 

8 

2,4 

1.5 

-63 

3 

3.2 

2,0 

-74 

18 

3  5 

2.5 

-81 

24 

3.6 

3.5 

-90 

35 

3.7 

4.5 

-95 

48 

3,8 

5  5 

-97 

60 

4.0 

7.5 

-99 

83 

4.1 

9.5 

-100 

100 

4.2 

11.5 

-100 

117 

4.2 

13.5 

-100 

130 

4.2 

15 

-100 

140 

4.3 

20 

-100 

135 

25 

-100 

128 

30 

-100 

110 

4.7 

oo 

-100 

14 

15 

The  4-meter  sheath  radius  can  be  estimated  by  equating  the  charge  in  the 
sheath  to  the  Laplacian  charge  on  the  probe; 
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Q  ~  4rtR^/3  Ne  ~  4jt  Eq  r  V 

which  gives  a  sheath  radius  of  3.7  meters.  The  current  can  be  estimated  as  the 
large  sheath  limit  of  the  Langmuir-Blodgett  cold  cathode  current  between  two 
concentric  spheres; 

I  =  0.8551  (4;teo)  (32e/81m)^^  {r/R)3  2 

which  gives  95  mA  for  the  parameters  of  this  problem.  As  this  is  far  greater  than 
the  plasma  thermal  current  to  the  sheath  surface  (1  mA),  the  sheath  expands  by 
plasma  erosion.  Equating  the  current  to  the  erosion  rate  gives  an  expansion 
rate  of  about  3  cm/psec.  Other  scaling  relations  can  also  be  defined.  We 
predict  a  time  of  about  0.01  seconds  to  expand  to  an  equilibrium  sheath. 


4.3  DYNAPAC  SIMULATION 

To  demonstrate  the  capability  of  performing  a  dynamic  particle-tracking  calculation, 
we  used  DynaPAC  to  calculate  in  three  dimensions  the  Gilbert  result  for  transient 
current  collection  by  a  1-foot  radius  sphere  biased  to  -100  kV.  The  procedure  was: 

1.  Initial  ions  were  placed  in  the  grid  to  mimic  a  uniform  density  of  10’^  m’^. 

2.  The  potential  was  initialized  at  a  time  appropriate  to  the  sphere  potential  at 
0.5  psec  using  barometric  electrons. 

3.  A  particle  pushing  timestep  of  0.5  psec  was  taken. 

4.  The  sphere  potential  was  advanced  0.5  psec,  and  one  major  potential 
iteration  performed.  The  -100  V  contour  is  taken  as  the  "sheath"  radius. 

5.  Return  to  3. 

Figure  4.1  shows  the  initial  potentials  with  -28  kV  on  the  sphere.  The  diamond¬ 
shaped  zero  contour  is  an  artifact;  contours  are  rounded  down  to  the  -10  V 
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contour  at  3.0  meters.  Analytically,  we  estimate  the  sheath  radius  to  be  a  bit 
over  2.4  meters,  which  is  the  location  of  the  -100  V  contour. 


Figure  4.2  shows  a  cut  of  the  ion  positions  after  the  first  0.5  psec  push.  Note 
that  a  uniform  ion  density  requires  a  slightly  nonuniform  ion  spacing. 

Figure  4.3  shows  the  second  set  of  potentials.  The  sphere  is  at  -49  kV,  and  the 
-10  V  contour  has  moved  to  3.5  meter  radius.  Figure  4.4  shows  the  second  set 
of  ion  positions,  indicating  visible  ion  movement  at  up  to  1 -meter  radius. 


Figure  4.5  shows  the  contours  at  3.5  psec.  The  10  V  contour  has  now  entered 
the  coarsest  (2-meter  resolution)  grid.  The  corresponding  ion  positions  are 
shown  in  figure  4.6.  Note  how  much  less  the  ions  from  the  corners  of  the  fine 
grid  have  moved  than  the  ions  from  the  nearest  grid  boundary.  The  sphere  is 
now  beginning  to  collect  ions  from  the  0.5  meter  resolution  grid. 
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•6  00E»03 
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.2  00E»03 
■1  OOE*03 
1  006*02 
■1  0OE*01 
OOOE.OO  (0) 


Figure  4.1  Potentials  around  sphere  at  0.5  psec. 
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Figure  4.2  Partfcle  positions  after  the  first  0.5  psec  push 
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Figure  4.3  Potentials  around  sphere  at  1 .0  psec. 


Figure  4.4  Particle  positions  after  1 .0  psec. 
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Figure  4.5  Potentials  around  sphere  at  3.5  psec. 


Figure  4.6  Particle  positions  after  3.5  jisec. 


Figures  4.7  through  4.10  show  potentials  and  particle  positions  at  9.5  and  15 
psec.  Ions  that  originated  just  inside  3  meters  are  just  entering  the  innermost 
grid  at  15  psec. 
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Figure  4.7  Potentials  around  sphere  at  9.5  psec. 


Figure  4.8  Particle  positions  after  9.5  psec. 
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Figure  4.9  Potentials  around  sphere  at  1 5  psec 
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Figure  4.10  Particle  positions  after  15  psec. 


Table  4.1  shows  the  DynaPAC  results  in  comparison  with  the  Gilbert  results. 
DynaPAC  shows  a  somewhat  faster  rise  in  current  to  a  somewhat  higher  peak 
value.  Also,  the  sheath  radius  is  somewhat  smaller.  This  is  consistent,  as  a 
smaller  sheath  leads  to  higher  electric  fields.  A  possible  reason  for  the 
difference  between  Gilbert  and  DynaPAC  is  that  the  DynaPAC  potentials  were 
not  fully  converged  at  each  timestep.  The  total  charge  collected  in  the 
DynaPAC  calculation  corresponds  to  the  charge  in  a  2.75  meter  radius 
spherical  region,  which  is  consistent  with  the  final  particle  positions  (figure 
4.10). 
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Table  4,3 


Time 

Potential 

Current 

Sheath  radius 

(M.sec) 

(kV) 

(mA) 

(m) 

0,5 

-28 

Gilbert 

DynaPac 

Gilbert 

DynaPAC 

4 

.2 

1.5 

2.4 

1.0 

-49 

8 

5.4 

2,4 

2.8 

1,5 

-63 

13 

15,4 

3.2 

3  2 

2.0 

-74 

18 

23.3 

3.5 

3.3 

2  5 

-81 

24 

31.8 

3.6 

3  4 

3.0 

-86 

40.9 

3.5 

-90 

35 

44.6 

3  7 

3  5 

4.0 

-93 

71.3 

4.5 

-95 

48 

87.3 

3  8 

3  6 

5.0 

-96 

38-8 

5.5 

-97 

60 

90.7 

4.0 

3  7 

6.0 

-98 

98.1 

6.5 

-99 

96.6 

7.0 

-99 

108.1 

7  5 

-99 

83 

96.1 

4.1 

3.8 

8.0 

-100 

133.7 

8.5 

■100 

122.7 

9.0 

-100 

72.1 

9  5 

-100 

100 

152.2 

4  2 

3.9 

10.0 

-100 

148.2 

10.5 

-100 

132.2 

11.0 

-100 

96.1 

11.5 

-100 

117 

156.2 

4.2 

4,0* ** 

12.0 

-100 

136.2 

12,5 

-100 

126.2 

13.0 

-100 

138.2 

13.5 

-100 

130 

120.2 

4.2 

4.1" 

14.0 

-100 

156.2 

14.5 

-100 

184.2 

15.0 

-100 

140 

68.1 

4.3 

4,2 

20.0 

-100 

135 

25,0 

-100 

128 

30.0 

-100 

110 

4.7 

Inf 

-100 

14 

15 

*  the  sheath  begins  to  move  to  grid  1 .  a  lower  resolution  grid. 

**  taking  sheath  contour  in  grid  2. 
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5.  APPLICATION  TO  SPEAR-II 


This  work  appeared  in  the  Interim  Report  (30  September  1990)  for  this 
contract.  The  predictions  were  presented  immediately  prior  to  the  SPEAR-II 

chamber  tests. 
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5.  APPLICATION  TO  SPEAR  II 


During  SPEAR  il  development,  the  NASCAP/LEO  code  was  used  to  calculate 
the  equilibrium  sheath  structure  and  currents  to  the  payload.  The  equilibrium 
sheath  was  relatively  large,  and  currents  were  focussed  on  the  high-voltage 
probe.  It  was  not  anticipated  that  the  currents  calculated  by  the  equilibrium 
code  would  have  any  noticeable  effect. 

Later,  a  "frozen  ion"  approximation  was  put  into  NASCAP/LEO.  This  showed 
that  the  transient  sheath  would  be  far  smaller  than  the  equilibrium  sheath. 
However,  it  was  not  known  what  the  implication  of  this  would  be  for  the  currents 
to  the  payload. 

Calculations  summarized  in  the  previous  section  indicated  that  transient 
currents  could  exceed  equilibrium  currents  by  an  order  of  magnitude.  Prior  to 
the  ground  test,  the  sphere  results  were  scaled  to  SPEAR  II,  using  the  sheath 
size  calculated  with  the  "frozen  ion"  approximation.  These  scaled  results 
suggested  that  the  current  to  the  high-voltage  probe  would  be  sufficient  to 
disturb  the  probe  reading.  The  ground  test  results  clearly  indicated  that  the 
high-voltage  probe,  which  functioned  well  in  vacuum,  was  unable  to  measure 
the  transformer  secondary  voltage  in  the  presence  of  the  laboratory  plasma. 
DynaPAC  calculations  currently  under  way  give  the  time-dependent  incident 
ion  currents  to  the  various  payload  components. 


5.1  SCALING  SPHERE  RESULTS  TO  SPEAR  li 

(These  results  were  presented  at  Plum  Brook  on  June  17.  1990.) 

Table  5.1  shows  scaling  of  the  transient  sheath  results  to  predict  peak  current 
collected  by  the  SPEAR  II  high-voltage  probe.  The  SPEAR  II  equilibrium  and 
transient  sheath  radii  were  inferred  from  calculations  using  the  static  plasma 
interactions  code  NASCAP/LEO  v/ith  the  default  and  "frozen  ion"  space-charge 
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nuia:  Tne  SPEAR  II  currents  were  estimated  from  tne  cbsenvation  tha! 

he  peak  cu"eru  occurred  at  twice  the  "transit  time"  and  cad  a  value  ot  cne- 
‘ourth  the  sneath  charge.  The  results  suggest  that  m  the  ia-ooratory  the  prooe 
will  coiiecit  C.1  amperes  of  incident  ion  current.  The  total  plasma  current  to  the 
orobe  wi!  be  about  one  ampere  due  to  secondary  electron  emission.  Since  the 
rominal  circ  ,it  current  through  the  probe  is  0.1  ampere,  this  will  cause  the 
probe  reading  to  be  in  error 


Table  5.1  Estimated  peak  currents  to  the  SPEAR  II  hic  -vonage  probe 


Sphere 

Flight 

Lab 

Density 

1  A  10'’ 

1  10” 

1  10'‘ 

Equilibrium  sheath  radius 

15  m 

2  3m 

t  8  m 

Equilibrium  current 

0  014  A 

0  0016  A 

0  010  A 

Transient  sneath  radius 

4  m 

1  3  m 

1  0  m 

Transit  time  {2R  mR'eV)'  ^ 

7  5  usee 

2  4  usee 

1  6  usee 

Sheath  charge 

4  3  uCOui 

0  15  ucoul 

0  7  ucoul 

Q.'4: 

0  14  A 

0  015  A 

0  1  A 

5.2  PLUM  BROOK  RESULTS 

The  SPEAR  II  ground  test  at  Plum  Brook  demonstrated  that  the  high-voltage 
components  could  be  conditioned  and  operated  both  in  vacuum  and  with 
plasma  present.  The  plasma  density  was  about  ■'O-^  m'^,  under  which  condition 
we  predicted  that  the  high-voltage  probe  would  collect  about  an  ampere  of 
plasma  current  (assuming  a  secondary  electron  yield  of  about  10)  and  that  its 
reading  would  be  perturbed. 

Figures  5.1  and  5.2  show  the  circuit  diagnostics,  without  and  with  plasma,  for 
clean  80  kV,  50  psec  pulses.  It  can  be  seen  that  all  diagnostics  are  the  same 
except  for  the  curve  labeled  "transformer  secondary  current,"  which  is  actually 
(proportional  to)  the  current  flowing  through  the  low-voltage  portion  of  the 
probe.  The  plasma  causes  the  probe  current  to  drop  to  about  a  quarter  of  its 
vacuum  value,  recovering  somewhat  toward  the  end  of  the  pulse. 

The  probe  consists  of  ten  0  1  MQ  resistors  in  series.  H  -100  kV  is  applied  to  the 
negative  end,  and  plasma  current  1  is  injected  into  the  node  one  resistor  away 
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Figure  5.1  SPEAR  II  voltage  probe  response  (upper  curve)  in  vacuum. 


Figure  5.2  SPEAR  II  voltage  probe  response  (lower  curve)  in  plasma. 
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from  the  negative  end.  the  current  through  the  remainder  of  the  probe  is  0,1(1 -I) 
amperes.  Thus  an  injected  piasma  current  of  a  sensible  fraction  of  an  ampere 
Will  cause  the  probe  reading  to  be  very  low.  as  is  seen  in  me  experimental 
results.  Had  the  plasma  density  been  an  order  of  magnitude  less,  as  might 
have  been  expected  in  flight,  the  probe  reading  would  have  been  less  severely 
affected,  perhaps  dropping  on  the  order  of  2C  percent. 
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6,  CALCULATION  OF  CURRENTS  TO  SPEAR  II 


This  work  appeared  in  the  Quarterly  Report  for  1  October  through  31 
December  1990.  It  was  subsequently  presented  at  the  High  Voltage  Workshop, 
Port  Hueneme,  CA,  19-21  March,  1991. 

Note  that  these  calculations  proceeded  at  a  rate  of  one  one-microsecond 
timestep  per  day  on  a  Sun  Microsystems  SPARCStation  I.  Subsequent 
DynaPAC  recoding  led  to  at  least  a  factor  of  three  speed  increase,  and  transfer 
to  an  SGI  Iris  Indigo  to  another  factor  of  three.  Were  we  to  repeat  the 
calculation,  we  should  be  able  to  do  several  timesteps  per  day. 
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6.  CALCULATION  OF  CURRENTS  TO  SPEAR  11 


DynaPAC  was  successfully  used  to  calculate  time-dependent  currents  to  the 
SPEAR  II  high-voltage  components  under  ground  test  conditions.  The  plasma 
was  specified  as  an  Ar+  plasma  of  density  10^2  m-3  j^e  pulse  shape,  intended 
to  represent  a  SPEAR  II  50  psec  pulse,  was  given  by 

V(t)  =  -1 .26  X  1 05  [1  -  exp(-t/3  x  1 0'®)]  exp(-t/50  x  1 0-6) 

This  voltage  function  has  a  peak  of  -100  kV  at  about  9  microseconds. 

The  calculation  was  carried  out  in  timesteps  of  one  microsecond,  with  each 
timestep  consisting  of  a  particle  pushing  operation  followed  by  a  potential 
solution.  On  a  SPARCstation  1 ,  the  calculation  proceeded  at  a  rate  of  about  one 
timestep  per  day,  with  about  20  percent  of  the  time  devoted  to  particle  pushing 
and  80  percent  to  potential  solving. 

Figure  6.1  shows  the  ion  current  incident  on  the  various  components.  The 
klystrode  battery-pack  and  the  transformer  (which  are  electrostatically  well- 
shielded  by  struts  and  bulkheads)  have  peak  currents  of  about  12  milliamperes, 
each  occurring  about  6  microseconds  into  the  pulse.  The  current  to  the  high- 
voltage  probe  (which  influences  a  larger  volume  of  plasma)  rises  to  a  peak 
exceeding  30  milliamperes  at  about  12  microseconds.  The  total  current  (which 
includes  current  to  the  high-voltage  leads)  peaks  around  55  milliamperes  at 
about  8  microseconds.  All  currents  fall  off  slowly  after  reaching  their  peak 
values.  For  comparison,  an  equilibrium  calculation  for  these  conditions  gives  a 
total  payload  current  of  about  6  milliamperes. 

Figure  6.2  shows  the  sheath  distance  as  a  function  of  time.  We  define  the  sheath 
distance  to  be  the  distance  from  the  high-voltage  end  of  the  probe  to  the  -100 
volt  contour.  During  the  first  5  microseconds,  as  the  potential  is  applied,  the 
sheath  grows  rapidly  to  its  “frozen  ion”  value  of  1.3  meters.  During  the 
remainder  of  the  calculation  it  grows  at  a  slow,  constant  rate  of  about  4 
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centimeters  per  microsecond.  This  growth  by  erosion  takes  place  even  though 
the  potential  is  decreasing  slightly  during  this  time. 

Figure  6.3  and  6.4  show  the  potential  contours  and  particle  postions  at  3 
microseconds.  This  is  past  the  peak  currents  to  the  klystrode  battery  and  the 
transformer,  and  it  can  be  seen  that  their  confined  regions  have  been  depleted 
of  ions.  Ions  can  be  seen  moving  toward  the  high-voltage  probe. 

Figure  6.5  and  6.6  show  the  potential  contours  and  particle  position  at  9 
microseconds.  Complex  flow  patterns  are  set  up  as  ions  travel  to  the  probe  from 
the  neighborhood  of  grounded  surfaces  such  as  the  top  bulkhead,  the  doors, 
and  the  plasma  accelerator. 

Figure  6.7  and  6.8  show  the  potential  contours  and  particle  postions  at  21 
microseconds.  Most  of  the  current  to  the  probe  now  consists  of  particles 
originating  a  meter  or  more  from  the  probe  in  grid  number  2.  However,  a  large 
number  of  macroparticies  originating  in  the  subdivided  regions  near  the  object 
are  still  alive.  Ions  originating  in  the  low-field  region  at  the  lower-left  corner  of 
grid  number  3  remain  unmoved. 
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Figure  6.3  Potential  contours  about  SPEAR  II  3  ^sec  from  the  beginning  of 
the  pulse. 
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Time  step  #  3  Total  tracking  time-  3.00E-06  sec. 
Particles  within  5,800  <-  Z  <»  6,000 
Mesh  Size  -  1  ,OOOOE+O0  meters 
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Figure  6.4  Ion  macroparticle  positions  at  3  psec. 


59 


X-AXIS 


Data  Name:  POTG  Time-  9 
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Figure  6.5  Potential  contours  about  SPEAR  II  9  |isec  trom  the  beginning  of 
the  pulse. 
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Time  step#  9  Total  tracking  time=  9.00E-06  sec 
Particles  wrthin  5.800  <=  2  <=  6.000 
Mesh  Size  =  1  OOOOE-t-CO  meters 


Figure  6.6  Ion  macroparticle  positions  at  9  psec. 
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X-AXIS 


Data  Name:  POTG  Time=  21 
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Figure  6.7  Potential  contours  about  SPEAR  II  21  ^isec  from  the  beginning  of 
the  pulse. 
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Y-AXIS 


Figure  6.8  Ion  macroparticle  positions  at  21  psec. 
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7.  SPEAR-II  PRODUCTS  REPORT 


This  chapter  is  the  Spring  1992  revision  of  S-Cubed’s  contribution  to  the 
SPEARTI  Products  Report.  The  final  extensive  rewriting  was  perlormed 
under  contract  F19628-9T  C-0187,  entitled  'Space  System-Environment 
Interactions  Investigation."  The  work  was  also  supported  by  contract 
DNA001-87-C-0091,  entitled  "Enhancing  Surface  Breakdown  Strength 
Components  in  Space." 

An  earlier  version  of  this  report  appeared  in  the  Quarterly  Report  for  1 
January  through  31  March  1991. 
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DEVELOPMENT  AND  APPLICATION  OF  PLASMA  INTERACTION  MODELS 
FOR  REALISTIC  HIGH  POWER  SPACECRAFT 


Introduction 

Recognizing  that  an  understanding  of  the  ambient  plasma  interactions  with  high  voltage,  pulsed 
power  systems  was  essential  to  extending  the  SPEAR  technology  to  other  systems  and 
environments,  two  existing  S-Cubed  plasma  simulation  codes,  NASCAP/LEO  and  Gilbert,  were 
challenged  to  make  SPEAR  predictions.  Since  neither  of  these  codes  could  perform  simulations 
which  were  both  dynamic  and  fully  three-dimensional,  the  SPEAR  program  contributed  to  the 
development  of  the  DynaPAC  {Dynamic  Plasma  Analysis  Code)  code*’' 

NASCAP/LEO  is  a  three-dimensional  computer  code  developed  by  the  NASA/Lewis  Research 
Center^'  for  the  study  of  high-voltage  plasma  interactions  and  the  design  of  high  voltage  systems 
in  LEO.  It  has  been  successfully  applied  'o  numerous  laboratory  and  spaceflight  experiments.  " " 
NASCAP/LEO  was  successfully  used  to  predict  plasma  currents  and  floating  potentials  for  SPEAR 
I,  whose  geometry,  while  truly  three-dimensional,  was  relatively  simple.  The  SPEAR  II  mock-up 
and  chamber  tests  showed,  for  the  first  time,  that  NASCAP/LEO  could  predict  sheath  plasma 
currents  for  a  complicated  geometry. 

Gilbert  is  a  general-purpose,  two-dimensional,  plasma  and  electrostatic  analysis  code  developed 
internally  by  S-Cubed.  The  Gilbert  calculations  performed  for  the  SPEAR  II  program  provided 
important  quantitative  assessments  of  the  validity  of  space  plasma  simulations  using  laboratory 
chambers.  Experiments  showed  that  the  computational  techniques  used  in  Gilbert  were  accurate 
enough  to  be  used  in  the  design  of  independent  components  for  space-based  power  systems. 
Pending  development  of  a  3-D  dynamic  code,  transient-to-equilibrium  current  scaling  relations 
developed  using  Gilbert  were  coupled  with  NASCAP/LEO  sheath  calculations  in  the  "equilibrium” 
and  "frozen-ion"  approximations  to  make  successful  pretest  estimates  of  enhanced  transient 
currents  to  the  SPEAR-I!  payload. 

DynaPAC  is  a  fully  three-dimensional,  dynamic  plasma  interaction  code,  which  includes,  through 
models  developed  and  validated  as  a  part  of  the  SPEAR  program,  much  of  the  knowledge  of 
plasma  physics  and  system/plasma  interactions  gained  during  the  SPEAR  program.  The  code 
allows  the  SPEAR  I  and  II  laboratory  data  to  be  extended  to  orbital  plasma  conditions  for  systems 
using  different  geometries. 

Applications 

NASCAP/LEO  models  high  voitage/piasma  interactions  in  three  dimensions  and  is  CAD/CAM 
compatible.  NASCAP/LEO  will  be  released  domestically  by  NASA/Lewis  Research  Center  in  the 
near  future.  When  released,  the  code  will  assist  researchers  in  the  study  of  high  voltage 
spacecraft/'plasma  interactions  and  spacecraft  designers  in  the  design  of  high  voltage  spacecraft 
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operating  in  LEO. 


Gilbert  is  a  two-dimensional  computer  code  for  the  modeling  of  static  and  dynamic  high 
voltage/plasma  interactions.  The  code  assists  researchers  in  the  study  of  high  voltage 
spacecrafVptasma  interactions  and  spacecraft  designers  in  the  design  of  high  voltage  spacecraft 
for  use  in  LEO.  Under  the  SPEAR  program,  the  code  was  applied  to  (1)  simulations  of  high 
voltage  breakdown  experiments  performed  at  Maxwell  Laboratories,  Inc..  (2)  modeling  of  the 
time-dependent  currents  to  the  SPEAR  II  payload  measured  in  a  smalt  space  simulation  chamber, 
and  (3)  modeling  of  the  sheath  expansion  around  and  the  transient  currents  to  a  high  voltage 
sphere. 

DynaPAC  allows  plasma  interaction  specialists  to  perform  realistic  analyses  with  direct  application 
to  engineering  problems.  Its  current  capability  is  illustrated  by  its  application  to  SPEAR  H.  In  the 
immediate  future,  DynaPAC  will  be  used  primarily  as  a  development  workbench  for  algorithms 
and  approximations  appropriate  to  different  regimes  of  spacecraft/plasma  interactions.  Eventually, 
a  version  of  DynaPAC.  containing  a  selection  of  validated  and  documented  plasma  formulations, 
will  be  released  for  engineering  use. 

In  addition  to  SPEAR  II,  many  other  space  experiments  and  space  power  system  analyses 
require  a  three-dimensional,  dynamic  treatment.  Some  examples  are  1)  spacecraft  charging  by 
an  electron  beam  at  an  arbitrary  angle  to  the  earth’s  magnetic  field,  2)  neutral  particle  beam 
charging,  3)  beam  plasma  interactions.  4)  beam/surface  interactions,  5)  ionization  instabilities 
within  electron  sheaths,  and  5)  ionization  breakdown  within  ion  sheaths  including  cathode  surface 
effects,  to  name  just  a  few. 

Prior  State-of-the-Art 

The  SPEAR-1  laboratory  and  flight  experiments,  along  with  three-dimensional  computer  modeling, 
demonstrated  the  capability  to  predict  steady  state  interactions  between  geometrically  simple  high 
voltage  spacecraft  and  the  space  plasma.  However.  (1)  the  ability  to  make  calculations  for 
spacecraft  which  are  truly  geometrically  complex  had  not  been  demonstrated,  and  (2)  the 
time-dependent  response  of  the  space  plasma  to  the  high  fields  and  voltages  associated  with 
pulsed  power  systems  such  as  SPEAR  II,  and  the  associated  dynamic  spacecraft  charging,  had 
not  been  fully  investigated.  Processes  not  adequately  modeled  include  formation  of  the  space 
charge  sheath,  current  flow  in  the  quasi-neutral  pre-sheath,  breakdown  phenomena,  plasma 
kinetics,  ionization  processes,  and  the  effect  of  dynamic  processes  on  spacecraft  charging. 

These  inadequacies  became  apparent  in  trying  to  make  plasma  interaction  predictions  for  the 
SPEAR  II  high  voltage  system.  Equilibrium  sheath  calculations  (e.g.,  using  NASCAP/LEO)  gave 
very  different  results  from  sheath  calculations  using  short  time  approximations.  The  plasma 
currents  to  the  high  voltage  components  could  therefore  not  be  calculated  with  confidence. 
Existing  Particle  In  Cell  (PIC)  codes  had  neither  the  generalized  geometry  necessary  for  realistic 
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systems  nor  the  sophisticated  algorithms  required  to  make  simulations  of  the  SPEAR  II  payload 
possible  in  a  reasonable  amount  of  computer  time.  Two  dimensional  dynamic  codes,  such  as 
Gilbert,  could  give  only  a  hint  of  the  detailed  information  needed  for  the  design  of  complex 
three-dimensional  systems. 

innovative  Concepts 

NASCAP/LEO  is  the  first  code  to  combine  plasma  sheath  physics  with  fully  CAO/CAM  compatible, 
three-dimensional  geometry.  Thus,  it  was  the  first  code  with  the  ability  to  model  a  payload  with 
the  geometrical  complexity  of  SPEAR  11. 

A  "frozen  ion"  capability  was  developed  and  incorporated  into  NASCAP/LEO  in  order  to  provide 
estimates  of  very  early  time  plasma  sheath  currents.  In  the  frozen  ion  approximation,  the  time 
scales  are  assumed  to  be  short  enough  that  inertia  prevents  the  ions  from  moving  under  the 
influence  of  the  fields.  This  was  a  stop-gap  method  for  calculating  time-dependent  currents  to  the 
SPEAR  II  payload  prior  to  the  development  of  the  DynaPAC  code. 

The  Gilbert  code  uses  grids  created  with  commercial  finite  element  programs  (such  as  PATRAN 
This  gives  the  ability  to  use  grids  with  variable  resolution  in  order  to  model  processes  in 
detail,  where  needed,  while  still  modeling  extended  range  interactions.  Higher-order  finite 
elements  with  continuous  electric  fields  and  a  third-order  particle  tracking  algorithm  are  used  in 
order  to  obtain  the  accuracy  needed  for  time-dependent  calculations.  The  code  allows  choosing 
the  appropriate  plasma  description  for  the  problem  under  consideration.  Models  based  on  particle 
tracking,  analytic  formulas,  and  hybrid  methods  are  available.  The  code  is  written  with  advanced 
programming  and  numerical  techniques  to  take  advantage  of  modern  computational  capabilities 
and  to  promote  ease  of  use,  ease  of  modification,  and  data  interchange  with  other  codes. 

In  order  to  model  SPEAR  laboratory  experiments  (described  later  in  this  chapter)  an  internal 
boundary  condition  to  mimic  the  presence  of  a  grid  (made  of  ordinary  window  screen)  was 
developed  and  incorporated  into  Gilbert.  It  was  shown  that  the  mean  potential  of  the  screen  was 
related  to  the  electric  field  discontinuity  across  the  screen. 

The  data  from  space  simulation  chamber  tests  of  the  SPEAR  II  payload  revealed  the  importance 
of  transient  plasma  phenomena  to  high  voltage  power  systems.  DynaPAC,  the  first 
three-dimensional,  dynamic,  plasma  interaction  code  which  resolves  the  complex  geometry  of 
realistic  systems,  was  developed  as  a  part  of  the  SPEAR  il  program.  This  code  was  written  with 
advanced  programming  and  numerical  techniques  to  take  advantage  of  modern  computational 
capabilities  and  to  promote  ease  of  use,  ease  of  modification,  and  data  interchange  with  other 
codes.  A  screen-handler  utility  provides  for  interactive  input  file  generation.  DynaPAC’s  DataBase 
Manager  is  a  programmer-friendly  way  of  allocating,  storing,  and  retrieving  large  blocks  of 
gridded  or  otherwise  structured  data.  DynaPAC  pioneers  the  use  of  higher  order  elements  which 
produce  strictly  continuous  electric  fields  and  potentials.  Arbitrarily  nested  grids  accommodate 
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simulations  of  complex  systems  with  extended-range  plasma  interactions.  Graphical  display  is 
available  using  a  variety  of  graphics  interfaces,  including  the  X-Windows  protocol,  which  allows 
display  from  a  remote  host. 

Description  of  the  NASCAP/LEO  Steady  State  Plasma  Sheath  Code 

NASCAP/LEO'^  *'is  a  three-dimensional,  finite-element  code  based  on  a  cubic  grid  structure.  For 
NASCAP/LEO,  a  spacecraft  is  defined  as  a  boundary  surface  element  representation  using  any 
standard  finite  element  preprocessor.  NASCAP/LEO  places  the  spacecraft  within  a  cubical  grid 
structure,  and  geometrical  descriptions  and  coupling  matrices  are  constructed  for  cubic  zones 
containing  the  object  surfaces.  Local  subdivision  of  the  basic  grid  is  used  to  resolve  critical 
regions,  and  nested  outer  grids  are  used  to  include  a  large  volume  of  space. 

NASCAP/LEO  was  used  to  calculate  the  plasma  current  distributions  to  the  SPEAR  II  payload. 
Such  calculations  are  based  on  the  concept  that  a  high  voltage  object  in  a  dense  plasma  forms 
a  sheath  within  which  the  plasma  is  highly  disturbed,  and  outside  of  which  the  plasma  is 
quiescent.  All  of  trie  computer  codes  described  in  this  chapter  solve  Poisson's  equation 

-CQ^<t>=p  (10‘1) 

in  various  ways,  where  $  is  the  potential,  to  is  the  permittivity  of  free  space,  and  p  is  the  space 
charge  density.  In  order  to  achieve  short  computational  times,  NASCAP/LEO  does  not  use 
particle  trajectory  information  to  calculate  the  space  charge  density  appearing  in  Eq.  (10-1). 
Instead,  a  nonlinear  analytic  expression  for  the  charge  density  as  a  function  of  the  local  potential 
and  electric  field,  based  upon  a  spherical  sheath,  is  used: 


p/eo=  -(4>M  d)(1  -  \m  !C((j),E))/(1 

C{<l>.E)-|0/4>|((RJr)2-1j  nO-2) 

(R^r)2=2.29|EAo/ej'“^|e/4)| 

p  =  space  charge  (coul-m  ’] 

Ca  =  8.854x10^^  [farad-m  ’l 
plasma  Debye  length  (m] 

6  =  plasma  temperature  [eV] 

0  =  local  space  potential  [volts] 

E  =  local  space  electric  field  [volts-m  ’] 

The  sheath  structure  is  usually  spherical  where  the  charge  density  most  greatly  influences  the 
potential,  so  Eq.  (10-2)  is  a  good  approximation  for  the  LEO  environment. 

Using  Eq.  (10-2)  for  the  space  charge  within  the  sheath  region,  together  with  the  appropriate 


68 


boundary  conditions  at  surfaces,  NASCAP/LEO  solves  Eq.  (10-1)  to  determine  the  electrostatic 
potential  about  a  spacecraft  having  a  highly  complex  geometry.  Knowing  the  potential  field,  the 
sheath  boundary  can  then  be  identified,  and  plasma  ions  and  electrons  crossing  the  sheath 
boundary  can  be  tracked  to  determine  the  current  distribution  on  the  spacecraft. 

Description  of  the  Gilbert  2-D  Dynamic  Plasma  Code 

Gilbert  is  a  two-dimensional  (R-Z  or  X-Y)  finite  element  computer  code.  It  solves  Poisson’s 
equation  (Eq.  10-1)  and  performs  particle  generation,  tracking,  and  deposition  on  a  grid  wnose 
elements  may  be  linear  triangles,  bilinear  quadrilaterals,  and/or  biquadratic  quadrilaterals.  A 
commercial  finite-element  preprocessor  (such  as  PATRAN  ‘’°‘)  is  used  to  generate  the  grid. 
Sequences  of  nodes  can  be  specified  for  use  in  assigning  electrostatic  boundary  conditions  or 
as  current  sources.  Elements  can  be  specified  as  empty-space  or  dielectric  filled  (with  optional 
conductivity). 

The  code  is  highly  modular  and  flexible.  Preprocessors  are  used  to  interpret  the  finite  element 
grid  input  and  to  generate  initial  particle  distributions.  The  main  analysis  routine  solves  problems 
either  by  time-stepping  or  by  iteration  on  a  nonlinear  problem.  Postprocessors  display  potential 
contour  plots,  particle  scatter  plots,  time  history  plots,  and  generate  and  display  particle 
trajectories. 

Description  of  the  DynaPAC  Three-Dimensional  Dynamic  Plasma  Code 

DynaPAC  (Dynamic  Plasma  Analysis  Code)  is  designed  to  perform  static  or  dynamic  plasma 
calculations  for  geometrically  complex  problems.  It  is  written  to  take  advantage  of  modern 
techniques  for  input  generation,  problem  solution,  and  visualization.  The  core  modules  of 
DynaPAC  allow  the  user  to: 

(1)  Define  the  spacecraft  geometry  and  the  structure  of  the  computational  space, 

(2)  Solve  the  electrostatic  potential  about  the  spacecraft,  with  flexible  boundary 
conditions  on  the  spacecraft  surfaces  and  a  space  charge  computed  either  fully 
by  particles,  fully  analytically,  or  in  a  hybrid  manner,  and 

(3)  Generate,  track,  and  otherwise  process  representative  macroparticles  of  various 
species  in  the  computational  space. 

The  core  modules  are  designed  to  have  the  maximum  practical  user  control  and  to  facilitate  the 
incorporation  of  new  or  modified  algorithms.  Preprocessors  are  provided  to  set  boundary 
conditions  and  generate  input  files  in  a  modern,  screen -oriented  way.  Similarly,  screen-oriented 
postprocessors  are  provided  for  graphical  and  textual  data  display. 
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Spacecraft  geometrical  definitions  are  done  using  standard  finite  element  preprocessors,  such  as 
PATRAN  Among  the  advantages  of  this  approach  are  that  the  geometry  can  be  realistically 
represented,  and  that  finite  element  models  of  a  spacecraft,  constructed  for  other  purposes,  can 
be  adapted  for  plasma  calculations.  The  computational  space  around  the  spacecraft  is 
constructed  interactively  using  the  GridTool  module.  Arbitrarily  nested  subdivision  allows 
resolution  of  important  object  features  while  including  a  large  amount  of  space  around  the 
spacecraft,  A  high-order,  finite  element  representation  for  the  electrostatic  potential  assures  that 
electric  fields  are  strictly  continuous  throughout  space.  The  electrostatic  potential  solver  uses  a 
conjugate  gradient  technique  to  solve  for  the  potentials  and  fields  on  the  spacecraft  surfaces  and 
throughout  the  surrounding  space.  Space  charge  options  presently  incorporated  include  Laplacian 
(appropriate  in  the  absence  of  plasma),  equilibrium  sheath  (appropriate  to  timescales  of 
milliseconds  or  longer),  "frozen  ions"  (appropriate  to  the  submicrosecond  stage  of  a  negative 
transient  pulse),  "mobile  ions  -  barometric  electrons"  (appropriate  to  the  SPEAR  11  case  of  several 
microsecond  time  scale  response  to  a  negative  pulse),  and  "full  PIC*  (appropriate  to  nanosecond 
timescales). 

Particle  tracking  is  used  to  study  sheath  currents,  to  study  particle  trajectories,  or  to  generate 
space  charge  evolution  for  dynamic  calculations.  Macroparticles  can  be  generated  at  either  a 
"sheath  boundary"  or  throughout  all  space.  Particles  are  tracked  for  a  specified  amount  of  time, 
with  the  time  step  automatically  subdivided  at  each  step  of  each  particle  to  maintain  accuracy. 
The  current  to  each  surface  cell  of  the  spacecraft  is  recorded  for  further  processing. 

The  strength  of  DynaPAC  lies  in  the  wide  range  of  length  scales,  time  scales,  and  physical 
phenomena  it  can  handle.  Arbitrarily  nested  subdivision  allows  the  modeling  of  small  geometrical 
features  within  large  scale  problems  (e.g.,  features  a  few  centimeters  in  size  in  sheaths  several 
meters  in  dimension).  By  suitable  choice  of  algorithms,  a  user  can  model  equilibrium  problems, 
problems  dynamic  on  an  ion  timescale  (microseconds),  or  problems  dynamic  on  an  electron 
timescale  (nanoseconds).  Trajectories  are  calculated  taking  account  of  user-specified  magnetic 
fields,  and  it  is  planned  to  introduce  a  procedure  to  obtain  potential  ana  space  charge  fields 
consistent  with  magnetic  effects. 

Gilbert  Application  to  High  Voltage  Breakdown  Experiments 

The  SPEAR  program  supported  a  series  of  experiments  at  Maxwell  Laboratories,  inc.,  to  study 
systematically  the  influence  of  a  background  plasma  in  promoting  the  breakdown  of  representative 
high  voltage  fixtures  ’'  Figure  10-1  shows  a  typical  test  fixture,  consisting  of  a  1  cm  diameter 
electrode  protruding  from  a  3  cm  diameter  insulator,  which  was  embedded  in  a  7  cm  diameter 
ground  plate.  In  vacuum,  this  fixture  would  hold  off  a  few  kilovolts  indefinitely  with  the  center 
electrode  at  a  negative  potential,  but  in  a  plasma  the  Insulator  surface  would  flash  over  in  a  short 
time. 
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Figure10-1 


Typical  axisymmetric  test  fixture  for  high  voltage  breakdown  experiments.  Two 
alternative  discharge-mitigation  screens  are  shown.  The  heavy  solid  line  indicates 
the  fixture  outline  as  seen  in  Figures  10-2  through  10-4. 
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if  the  plasma-induced  discharge  resulted  from  electric  field  modification  by  u  . 

collected  on  the  dielectric  surface,  we  would  expect  the  hold-off  time  to  be  inversely  prGp:.;f  or  a 
to  the  plasma  current.  Supposing  this  to  be  the  case,  a  hemisphere  of  ordinary  a.'sm  f'sc;  ,*.1;  co,- 
screen  was  mounted  on  the  fixture  in  order  to  reduce  the  total  ion  current  coilectea  Ca:  .. 

showed  that  the  screen  would  reduce  the  ion  current  by  approximately  a  factor  ot  five 
Experimentally,  the  screen  was  found  to  increase  the  holdoff  time  from  30  ms  to  30  s  at  -2  5  s 
and  from  15  ms  to  5  s  at  -5.0  kV  in  a  background  plasma  density  of  4xi0^  cm  ^  Tnis  incic-aso  c 
the  holdoff  time  was  too  great  to  be  explained  solely  in  terms  of  the  total  plasma  current  to  r- 
fixture. 

Calculations  of  the  plasma  ion  trajectories,  performed  with  the  Gilbert  code,  were  used  to  ai.e.r 
to  understand  this  phenomenon.  Figure  10-2  shows  the  trajectories  of  ions  impinging  cn  ;no 
original  fixture  (without  the  screen).  It  is  apparent  that  ionu  impact  the  entire  surfaca  of  tn.c  j  a 
including  the  sensitive  "triple  point"  area  where  the  insulator  contacts  the  high  voltage  eiectroou 

In  order  to  study  the  ion  trajectories  in  the  presence  of  the  screen,  it  was  necessary  !c  ou'  e  cp 
an  internal  boundary  condition  that  would  mimic  the  presence  of  the  screen.  It  was  found  thai  tne 
mean  potential  of  the  screen  was  related  to  the  electric  field  discontinuity  across  the  scree  1 
Figure  10-3  shows  the  ion  trajectories  with  the  "normal”  screen  present.  The  ion  current 
focused  onto  the  top  of  the  center  electrode,  with  no  calculable  ion  current  to  the  triple  point 
Thus,  the  effect  of  the  screen  was  to  steer  plasma  ions  away  from  the  triple  point,  m  add'tion  to 
reducing  the  total  ion  current  to  the  fixture. 

In  order  to  confirm  this  finding,  a  further  calculation  was  performed  which  showed  tnat  if  tm: 
screen  was  configured  horizontally  at  its  attachment  point,  then  there  would  be  calculable  current 
to  the  triple  point  (Figure  10-4).  The  fixture  with  the  horizontal  screen  broke  down  m  40  rns  witr, 
an  applied  voltage  of  -15  kV  and  a  plasma  density  of  8x10^  cm^,  compared  with  750  nis  for  tire 
norma!  screen.  Additional  calculations  and  experiments,  involving  larger  fixtures  ano  vano.rc 
screen  shapes,  showed  good  correlation  between  rapid  breakdown  and  the  calculated  lon  curren; 
to  the  triple  point.  (See  Tablel0-1.)  These  studies  also  illustrated  that  trying  to  impfov-.; 
Dreakdown  characteristics  by  the  obvious  step  0^  enlarging  the  insulator  may.  in  fact.  p-Gmo'.,t 
brea.Kdown  by  reducing  the  electrostatic-inertial  insulation  of  the  triple  point.  These  results  she  ,' 
the  ut'iity  of  accurate  two-dimensional  potential  and  trajectory  calcjiations  m  tiesignir,; 
axisyrnmetne  high  voltage  fixtures  for  maximum  holdoff  time. 

Gilbert  Small  Chamber  Model  for  SPEAR  II 

The  SPEAR  II  payload  was  first  tested  with  a  background  plasma  in  a  small  vacuum  chamber  .51 
Maxwell  Laboratories,  Inc.  It  was  observed  that  the  10  ps  pulses  were  more  likely  to  cause 
breakdown  than  the  3  or  50  ps  pulses.  This  finding  motivated  a  computer  simulation  using  Giiben 
of  the  plasma  interactions  in  the  small  chamber  to  determine  if  plasma  dynamic  effects  might  me 
responsible  tor  this  phenomena. 
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Table  10-1.  Holdoft  time  for  various  insulator  widths  and  screen  shapes 


Insulator  Width 

Applied  Voltage 

Time  to  Breakdown 

Trajectories  to  Triple  Point 

1.0  cm  (n) 

-15  kV 

750  ms 

No 

1.0  cm  (h) 

-15  kV 

40  ms 

Yes 

1.0  cm  (n) 

-50  kV 

20  ms 

No 

2.0  cm  (n) 

-50  kV 

70  ms 

No 

4.5  cm  (1) 

-50  kV 

7  ms 

Yes 

4.5  cm  (s) 

-50  kV 

3  ms 

Yes 

Time  to  breakdown  tor  fixtures  with  differing  insulator  widths.  The  plasma  density  "'as 
8x10'  cm  '.  The  screen  shapes  are  (n)  normal  screen  (Figure  10-1),  (h)  horizontal  screen 
(Figure  10-1),  (1)  "long"  screen  (not  shown),  and  (s)  “short"  screen  (not  shown). 


Figure  10-5  shows  an  artist's  concept  of  the  two-dimensional  axisymmetric  representation  of  lee 
payload  used  in  the  calculation.  A  grid  was  constructed  to  represent  the  space  between  the  mocei 
and  chamber  walls  and  ion  macroparticles  were  placed  in  the  grid  to  represent  ar.  initial  uniform 
Ar  plasma  with  a  density  of  1x10'  cm  '.  The  applied  voltage  on  the  model  was  increased  with  .a 
risetime  appropriate  to  the  pulse  width  under  study,  and  the  ions  were  allowed  to  move  in  me 
resulting  fields.  Figure  10-6  shows  the  ion  macroparticle  positions  4  ps  into  the  puise.  The  ior';s 
are  seen  to  be  focused  on  the  high  voltage  end  of  the  voltage  divider,  precursive  to  the  resu'tr, 
later  obtained  in  three  dimensions  with  DynaPAC.  By  this  time,  about  half  the  plasma  ions  (whici  i 
initially  filled  the  chamber  with  a  nearly  uniform  dot  pattern)  are  gone.  The  chamber  will  be  almost 
completely  drained  of  plasma  ions  within  about  10  |xs.  (The  plasma  source  is  not  able  to  replenrsn 
the  plasma  on  so  short  a  timescale.) 

Figure  10-7  shows  the  calculated  current  to  the  payload  for  two  cases:  (1)  a  risetime  of  1.5  ps 
(characterizing  the  3  ns  and  10  ns  SPEAR  il  high  voltage  pulses)  and  (2)  a  risetime  of  3  ns 
(characterizing  the  50  ns  high  voltage  pulse).  The  peak  current  occurred  at  about  3  ns  into  the 
pulse  for  both  cases,  but  is  noticeably  higher  for  the  1.5  ns  risetime  pulse.  There  is  no  current 
beyond  10  ns,  as  the  chamber  is  drained  of  plasma.  The  suggested  interpretation  is  that: 

(1)  The  3  ns  pulse  does  not  lead  to  breakdown  because  the  pulse  is  over  when  ti  i  ,. 
current  peak  occurs. 

(2)  The  10  ns  pulse  is  more  likely  to  cause  breakdown  than  the  50  ns  pulse  f: 
the  higher  peak  incident  ion  current. 

(3)  Extending  the  pulse  beyond  10  ns  does  not  lead  to  breakdown  because  there  i$ 
no  further  ion  current  to  the  payload. 

The  overall  conclusion  from  these  calculations  is  that  a  small  chamber  may  produce  resuns 
unrepresentative  of  the  behavior  of  a  system  in  a  very  large  chamber  or  m  the  ionosphere. 
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Figure  10-5.  Artist's  representation  of  the  axisymmetric  mooel  of  SPEAR  II  used  in  the  two 
dimensional  plasma  dynamic  studies.  The  labels  indicate  the  peak  potentials 
applied  to  each  component:  -100  kV  on  the  probe  lead,  a  linear  gradient  on  the 
high  voltage  probe  from  -100  kV  at  the  probe  lead  to  0  at  the  front  plate,  and  on 
the  klystrode  section  (representing  a  cavity  partially  enclosed  by  several  struts) 
linear  gradients  from  0  at  either  end  to  -50  kV  at  the  center. 
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Figure  10  6.  Gilbert-calculated  ion  macropartide  positions  4  ps  into  a  SPEAR  n  l 
pulse. 
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Figure  10-7, 


3.0  ^  s  PULSE  RISE  TIME  -50p.s  PULSE  WIDTH 


Oe+0  2g-6  4G“6  OoG  8o*6  ie-5 


The  calculated  ion  currents  to  the  axisymmetric  model  for  1 .5  and  3  rise  time 
pulses  (see  text  for  explanation). 
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SPEAR  II  Geometrical  Model  for  NASCAP/LEO  and  DynaPAC 


Figure  10-8  shows  the  geometrical  mode!  of  the  SPEAR  I!  payload  used  ‘or  U'.e  f-.-ASCAt-'  .  ■  ;  ■ 
and  DynaPAC  calculations.  (The  identical  model  was  used  for  the  SPEAF^  II  mock-up  ;  A  ■> 
geometry  was  based  on  the  EUCLID  (an  engineering  araftmg  and  modeling  product  ■' i  rnc  ju 
constructed  by  Westinghouse  and  was  converted  to  a  PATRAN  (a  f;r„te  element  rncau  '  < 
code''°)  model  using  IGES  (Initial  Graphics  Exchange  Specification  protocol  ')  files  as 
intermediate  transfer  mechanism.  The  PATRAN  "Neutral  File"  (specifying  node  locations  -.e  j 
surface  element  connectivities,  materia!  numbers,  and  conductor  numbers)  is  read  oy  the  or; 
definition  interface  module  of  either  code,  which  places  the  object  in  the  grid  system  and  perrerms 
other  appropriate  preprocessing. 

Figure  10-9  shows  the  DynaPAC  grid  for  SPEAR  11,  illustrating  the  arbitrarily  nested  subdiv  s  i.,; 
capability.  The  grid  spacing  ranges  from  1  meter  (outermost  grid,  not  shown)  to  3  ern  in  the 
payload  interior.  The  NASCAP,  LEO  grid  had  a  basic  mesh  size  at  16  cm.  wit.n  subdivided  regie 'i.y 
down  to  4  cm,  and  surrounding  outer  grids  up  to  64  cm. 

NASCAP/LEO  Calculations  for  SPEAR  It 

Prior  to  the  development  of  DynaPAC,  NASCAP'LEO  was  the  only  plasma  interaction  c;;:n 
capable  of  predicting  sheath  currents  and  potentials  for  a  payload  with  the  geometrical  compic  ^n  . 
of  SPEAR  II,  It  was  recognized  that  the  "equilibrium"  plasma  treatment  used  in  f^ASCAP/LEO  .'.cr 
not  appropriate  to  the  SPEAR  II  high  voltage  pulse  lengths.  Nonetheless,  predictions  made  us  ' 
an  equilibrium  code  were  better  than  no  predictions  at  ail.  Figure  10-10  snows  the  NASCAP  i_b 
calculated  equilibrium  sheath  contours  tor  SPEAR  II. 

As  part  of  the  SPEAR  program,  a  "frozen  ion”  approximation  was  added  to  t-iASvCAP'LEO  In  t  - 
"frozen  ion"  approximation,  it  is  assumed  that  electrons  instantaneously  attain  an  ^jputiibn./. 
distribution,  but  ions  do  not  move.  The  space  charge  function,  eq.  {10-2).  is  thus  replaced  oy 


p/^o  =  0  4>^0 

p/eo  =  (0/4)(1  (J)^0 

This  approximation  (appropriate  up  to  about  1  ps  for  SPEAR  II  parameters)  prruvided  a  mca-  : 
for  calculating  sheath  sizes  during  the  short  SPEAR  il  high  voltage  pulses.  i.See  Figure  10  ■  r  . 
However,  it  did  not  provide  any  straightforward  means  of  calculating  the  transient  currents  to  i-'ie 
payload  or  the  time  dependence  of  the  currents  and  sheath  structure. 

The  first  question  to  be  addressed  using  NASCAP/LEO  was  the  floating  potential  of  the  payi'i.i-d 
during  high  voltage  operation.  It  was  not  sufficient  simply  to  assume  that  the  floatmc  pote  i!  ,T 
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Figure  10-9.  DynaPAC  grid  for  SPEAR  !!  (The  outmost  grid,  with  1  meter  mesh  spacing,  is 
omitted.) 


82 


Figure  10-10,  NASCAP/LEO  calculated  equilibrium  sheath  ("H"  contour)  about  SPEAR  II  under 
laboratory  conditions  {1x10®cm^;  1  eV).  Space  charge  densities  are  calculated 
using  formulas  appropriate  to  a  steady-state  plasma  sheath  (e.g.,  Langmuir- 
Biodgett).  The  shaded  area  represents  the  projection  of  the  payload  on  the 
contour  plans,  and  equipoential  lines  with  10  kV  spacings  appear  near  the 
payload.  The  contour  at  -100  volts  (marked  "H")  is  considered  the  sheath  contour 
because  the  potential  should  drop  to  zero  within  0.3  meters  of  the  -100  volt 
equipotential  surface.  Also  shown  are  contour  lines  at  -10  volts  (marked  ’T")  and 
at  zero  volts  (marked  "0").  The  payload  ground  surface  is  at  +100  volts.  The 
horizontal  and  vertical  axes  are  mai'ked  in  units  of  16.3  cm. 
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Figure  10-11.  NASCAP/LEO  calculated  "frozen  ion"  sheath  ("H"  contour)  about  SPEAR  II  under 
laboratory  conditions.  Space  charge  densities  are  calculated  assuming  thai  iun 
motion  is  negligible  (appropriate  to  times  up  to  a  few  microseconds)  Othecvvise, 
the  description  of  Figure  10-10  applies.  Note  that  the  transient  sheath  is  tar 
smaller  than  the  equilibrium  sheath,  leading  to  considerably  elevated  currents  m 
the  transient  regime. 
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calculating  the  electron  and  ion  currents  to  the  payload  for  several  assumed  values  of  floating 
potential,  NASCAP/LEO  predicted  that  the  equilibrium  floating  potential  would  be  less  than  1 50  V. 
This  prediction  was  confirmed  insofar  as  the  SPEAR  il  mockup  instrumentation  was  unable  to 
show  any  evidence  of  a  positive  steady  state  floating  potential. 

The  code  was  next  used  to  study  the  current  distribution  to  the  various  components  of  the  SPEAR 
II  payload.  For  the  mockup  test  under  positive  bias,  NASCAP/LEO  predicted  a  total  current  of 
-1.2  amperes  for  a  plasma  of  density  1x10®  cm  ^  and  temperature  1  eV.  The  distribution  of  this 
current  to  the  SPEAR  il  components  is  shown  in  Table  10-2. 


Table  10-2  Currents  to  SPEAR  il  Components 

SPEAR  II  Component 

Calculated  Current 

100  kv  portion  of  probe 

500  ma 

75  kv  portion  of  probe 

410  ma 

Klystrode  battery  canister 

87  ma 

50  kV  portion  of  probe 

75  ma 

Various  grounded  surfaces 

59  ma 

High  voltage  lead 

49  ma 

Transformer  secondary 

34  ma 

25  kV  portion  of  probe 

5  ma 

Klystrode  bushing 

0,2  ma 

NASCAP/LEO  calculation  of  currents  to  components  of  the  SPEAR  II 
mockup  under  positive  bias  conditions. 


The  NASCAP'LEO  calculations  predicted  that  nearly  all  of  the  current  would  be  incident  to  the 
high  potential  end  of  the  high  voltage  probe.  This  prediction  was  qualitatively  confirmed  by  visual 
observations  during  the  mockup  test  of  a  strong  optica!  glow  from  the  electron  bombarded  region. 
However,  the  question  of  how  dynamic  effects  would  modify  this  current  distribution  in  the 
negative  bias  case  remained  open. 

Transient  Current  Estimates  Using  Gilbert  and  NASCAP/LEO 

it  was  recognized  that  equilibrium  current  calculations  for  the  SPEAR  II  payload  could  not  be 
trusted  for  the  short  pulse  durations.  Without  a  three-dimensional  modeling  capability  it  was  not 
known  to  what  degree  the  transient  current  would  differ  from  the  equilibrium  current,  or  how  much 
time  was  required  to  achieve  the  equilibrium  configuration.  To  shed  some  light  on  these  topics, 
a  two-step  strategy  was  used.  First,  the  Gilbert  code  was  used  to  calculate  the  time-dependent 
current  to  a  sphere  with  parameters  similar  to  the  SPEAR  11  system.  From  the  results,  scaling 
relations  were  developed  relating  the  "frozen  ion"  sheath  size  to  the  peak  transient  current. 
Second,  NASCAP/LEO  was  run  to  calculate  the  size  of  the  "frozen  ion"  and  equilibrium  sheaths 
about  the  high  voltage  probe,  and  the  scaling  relations  were  applied  to  determine  the  peak 
transient  current. 
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The  sample  problem  consisted  of  a  0.3  m  radius  sphere  in  an  O*  plasma  wun  a  density  c!  ;  , 
cm'^.  The  potential  was  raised  to  -100  kV  with  a  1.5  ^is  risetime.  Using  spherical  proPe  itiojr.. 
the  equilibrium  sheath  around  the  sphere  was  calculated  to  have  a  radius  of  15  m,  ana  tne 
steady-state  ion  current  to  the  sphere  was  14  mA. 

The  Gilbert  calculation  was  run  for  30  us.  it  was  found  that  the  transient  sheatn  expands  rap  a  y 
to  a  radius  of  4.0  m  as  the  voltage  is  applied,  and  thereafter  expands  by  plasma  erosion  a;  me 
slow  rate  of  3  cm/ps,  reaching  a  radius  of  4.7  m  at  the  end  of  the  calculation.  The  current  to  toe- 
sphere  (Figure  10-12)  reaches  a  peak  of  140  mA  at  15  ps  into  the  pulse,  and  fans  off  sro  viy 
thereafter.  (These  calculations  gave  the  first  indication  that  the  SPEAR  II  incident  ion  currents 
might  be  an  order  of  magnitude  above  the  values  obtained  from  the  equilibrium  treatment.)  it  was 
estimated  that  10  ms  would  be  required  for  the  current  to  relax  to  the  Langmuir-Blcdgetr  ^ 
equilibrium  current  of  14  mA.  (Note  tnat  the  collected  current  correlates  inversely  with  the  sn^-pt.  i 
size.  This  is  analogous  to  the  space-charge-limited  current  collected  by  an  inner  sphere  fror.-  a 
concentric  outer  sphere.) 

To  develop  scaling  a  relation,  we  postulate  that  the  peak  current  is  given  py 

where  Q  is  a  characteristic  charge  and  x  a  characteristic  time.  The  characteristic  cha',,'; 
logically  taken  as  the  space  charge  in  the  sneath  ywhich  must  balance  the  surface  chatge  a 
probe); 

Q  M4r-/3)  ne  r|,  n  o  ¬ 

where  we  have  neglected  the  probe  volume  relative  to  the  sheath  volume.  The  characteristic  t;r 
may  be  taken  as  the  time  for  an  ion  to  travel  the  sheath  radius  m  a  uniform  field  given  by 
applied  voltage  over  the  sheath  radius: 

x=(2mR2/eV)^ 

For  the  Gilbert  calculation,  we  have  Q  =  4.3x10’''  coulombs,  and  x  =  7.3x10^  seconas,  so  tna:  t 
scaling  constant  takes  the  value 

a  =  'peak  '/Q  = 


and  the  scaling  relation  is 


Current  {milliamperesj 


]Table10-3  shows  "frozen  ion”  estimates  of  the  sheath  size  (calculated  by  NASCAP/LEO  for  tiio 
three  din iensional  representation  of  the  SPEAR  li  payload  under  space  and  laboratory  conditions), 
and  the  peak  transient  currents  to  the  high  voltage  probe  (predicted  by  scaling  the  Gilbert  results). 
To  obtain  the  actual  parasitic  current  in  the  circuit,  the  incident  ion  current  must  be  enhanced  by 
a  large  factor  to  take  into  account  the  secondary  electron  yield  of  the  incident  ions.  While  there 
is  a  great  deal  of  uncertainty  in  the  secondary  yield  value,  a  yield  of  25  would  be  a  reasonable 
estimateJ’^^This  gives  a  parasitic  current  exceeding  two  amperes  under  laboratory  conditions.  For 
comparison,  the  fully  dynamic  results  (obtained  from  DynaPAC,  as  described  below)  predicted 
a  total  peak  incident  ion  current  of  50  mA,  with  35  mA  on  the  probe,  for  a  parasitic  current  of  0,9 
amperes. 


Table  10-3.  Transient  Current  Estimates 

Gilbert 

Space 

Laboratory 

Density  [cm  ’j 

1  X  10^ 

1  X  10^ 

1  X  10® 

Equilibrium  Sheath  Radius  [m] 

15 

2.3 

1.8 

Equilibrium  Current  [mA] 

14 

1.6 

10 

Transient  Sheath  Radius,  Rg  [m] 

4.0 

1.3  m 

1.0  m 

Transit  Time,  t  [ps] 

7.3 

2.4 

1.6 

Sheath  Charge  [pcoul] 

4.3 

0.15 

0.7 

Q/4t  [mA] 

140 

15 

11 

Secondary  Electron  Yield 

25 

25 

25 

Peak  Pa.'asitic  Current  [A] 

3.6 

0.39 

2  6 

Estimates  of  peak  parasitic  current  (1)  calculated  by  Gilbert  for  a  0  3 
meter  radius  sphere;  (2)  estimated  for  space  conditions  (plasma  density  10‘ 
cm^)  using  the  NASCAP/LEO  "frozen  ion"  approximation;  and  (3)  similarly  estimated  tor 
laboratory  conditions  (plasma  density  10‘  cm’). 

The  predicted  current  levels  would  not  impair  the  high  voltage  system  operation,  but  v 
seriously  impact  the  ability  of  the  high  voltage  probe  to  measure  the  transformer  seconcary 
voltage.  Figure  10-13  shows  a  circuit  model  of  the  probe,  in  which  parasitic  plasma  current  is 
injected  one-tenth  of  the  way  from  the  high  voltage  end.  and  the  voltage  measurement  is  taken 
one-tenth  of  the  way  from  the  ground  end.  For  the  case  of  -100  kV  applied  to  the  probe,  this 
circuit  model  predicts  that  a  parasitic  current  of  0.5  amperes  will  result  in  the  measurement  being 
low  by  a  factor  of  two. 

A  comparison  of  the  probe  traces  obtained  from  the  space  simulation  chamber  tests  under 
vacuum  and  plasma  conditions  dramatically  illustrates  this  effect.  Figure  10-14  shows  diagnostic 
signal  traces  for  a  50  us,  80  kV  pulse  in  vacuum.  The  high  voltage  probe  trace  is  the  dark,  upper 
trace;  the  other  traces  are  diagnostics  taken  at  other  components  of  the  high  voltage  circuit  In 
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Figure  10-13.  Circuit  model  of  high  voltage  probe,  used  to  estimate  the  effect  of  parasitic  plasma 
current.  The  transformer  secondary  voltage  is  measured  as  lOxV^^,  =  (1 .  -1pia5„,a) 


Figure  10-1 5  the  trace  is  shown  tor  the  same  pulse  parameters  in  the  presence  of  a  10®  cm  ^  Ar 
plasma.  The  high  voltage  probe  reading  is  diminished  by  about  a  factor  of  two,  but  all  of  the  other 
diagnostic  signals  are  ind^Tinguishabie  from  the  vacuum  shot.  Our  conclusion  is  that  the  high 
voltage  circuit  is  operating  properly  in  the  presence  of  plasma,  but  the  high  voltage  probe  is 
reading  incorrectly. 

DynaPAC  Calculation  of  Transient  Currents  to  SPEAR  11  High  Voltage  Components 

All  the  calculations  described  above  either  use  an  equilibrium  or  "frozen  ion"  space  charge 
representation,  neither  of  which  adequately  model  the  dynamics  of  the  SPEAR  II  pulse,  or  else 
model  the  geometry  as  axisymmetric,  which  is  not  an  adequate  approximation  to  the  partially- 
enclosed  structure  of  SPEAR  II.  Estimates  of  the  peak  incident  ion  current,  the  division  of  this 
current  among  the  SPEAR  II  high  voltage  components,  and  the  timescale  of  the  current  obtained 
by  scaling  the  results  of  equilibrium  codes  based  on  results  for  poorly  representative  geometry 
required  confirmation  by  a  fully  3-D  dynamic  code.  For  this  reason  the  DynaPAC  code,  already 
under  development  for  Geophysics  Laboratory  (now  Phillips  Laboratory)  was  supported  by  the 
SPEAR  program. 

With  the  development  of  DynaPAC  for  the  SPEAR  fl  program,  it  was  possible  to  calculate  the  time 
dependent  ion  currents  to  the  various  SPEAR  II  high  voltage  components.  The  geometrical  model 
of  the  SPEAR  II  payload  was  the  same  as  was  used  for  the  NASCAP/LEO  calculations  described 
above.  Initially,  each  DynaPAC  grid  was  filled  with  a  regular  array  of  ion  macroparticles, 
representing  an  Ar  plasma  with  a  density  of  1x10®  cm'®.  At  each  timestep,  a  time-dependent 
voltage,  representing  a  100  kV,  50  ixs  pulse  with  3  ps  risetime,  was  applied  to  the  high  voltage 
components,  the  electrostatic  potential  field  was  calculated  using  the  known  ion  distribution  and 
an  analytic  expression  for  the  electron  charge  density,  and  the  ion  macroparticles  were  tracked 
in  the  new  field  for  the  specified  timeutep  length.  After  3  ps,  the  ion  macroparticles  assumed  the 
configuration  represented  in  Figure  10-16.  The  figure  shows  ion  voids  near  the  klystrode  battery 
pack  and  the  transformer.  This  indicates  that  substantial  ion  motion  has  taken  place  in  these 
regions,  which  had  high  fields  due  to  the  close  proximity  of  grounded  struts  and  bulkheads.  By 
contrast,  ions  at  large  distances  were  just  beginning  to  be  accelerated  toward  the  high  voltage 
probe,  as  indicated  by  the  curvature  of  the  originally  straight  lines  of  macroparticles. 

Figurel  0-1 7  shows  the  calculated  ion  currents  incident  on  various  SPEAR  II  payload  components. 
The  klystrode  battery  pack  and  transformer  (which  are  well  shielded  electrostatically  by  struts  and 
bulkhead)  had  peak  currents  of  about  12  mA  each,  which  occurred  approximately  6  ps  into  the 
pulse.  The  current  to  the  high  voltage  probe  (which  influences  a  larger  volume  of  plasma)  rose 
to  a  peak  exceeding  30  mA  at  about  12  ps.  The  total  incident  ion  current  (which  includes  currents 
to  the  high  voltage  leads)  peaked  at  around  55  mA  approximately  8  ps  into  the  pulse.  Ail  of  the 
currents  fell  off  slowly  after  reaching  their  peak  values.  (For  comparison,  NASCAP/LEO 
equilibrium  calculations  for  these  conditions  gave  a  total  payload  current  of  6  mA.  See 
Table  10-4.) 
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Figure  10-15.  Circuit  diagnostics  for  an  80  kV,  50  microsecond  SPEAR  II  pulse  with  plasma 
sources  operating.  The  transformer  secondary  voltage  measurement  is  about  naif 
the  value  measured  in  vacuum,  and  now/  appears  as  the  lower  curve.  Remaining 
circuit  diagnostics  are  unchanged.  DynaPAC  calculations,  anticipated  by 
NASCAP/LEO  and  Gilbert  results,  showed  that  the  voltage  divider  probe 
measurement  was  incorrect  in  the  presence  of  plasma  due  to  high  incident  lon 
currents  to  the  probe. 
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Figure  10-16.  Calculated  ion  macroparticle  positions  3  into  a  SPEAR  il  high  voltage  pulse 
(scale  marked  in  units  of  meters.) 
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Figure  10-17.  Calculated  incident  ion  currents  to  the  SPEAR  I!  high  voltage  probe,  klysirode 
battery  canister,  and  pulse  transformers. 
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Table  l0-4  Predicted  Incident  Ion  Currents  to  SPEAR  I! 

Equilibrium  Calculation 

6  ma 

Pre-Test  Prediction 

100  ma 

DynaPAC  Calculatkjn 

55  ma 

Taking  into  account  the  secondary  electron  yield,  the  predicted  parasitic  current  through  the  upper 
part  of  the  probe  is  nearly  one  ampere.  As  discussed  in  the  previous  section,  this  is  easily  enough 
current  to  cause  the  observed  erroneous  measurements  by  the  high  voltage  probe. 

A  rough  estimate  of  the  time  needed  to  approach  an  equilibrium  sheath  may  be  obtained  as  the 
time  required  to  drain  the  equilibrium  sheath  volume  (approximately  a  2  1  meter  radius  sphere) 
of  charge  at  a  mean  current  of  about  40  mA.  This  approach  gives  a  sheath  formation  time  of  1 50 
ps,  which  is  three  times  the  longest  SPEAR  II  high  voltage  pulse.  Therefore,  the  equilibrium 
formulation  is  never  appropriate  for  the  SPEAR  II  operating  conditions. 

Summary  of  Technical  Achievements  and  Capabilities 

This  chapter  has  discussed  the  development  of  three  plasma  interaction  model  computer  codes 
and  their  application  to  the  SPEAR  program.  The  two  existing  codes,  NASCAP/LEO  and  Gilbert, 
were  both  enhanced  under  the  SPEAR  program,  and  their  ability  to  make  useful  and  non-obvious 
predictions  was  dem,onstrated.  The  new  code.  DynaPAC,  demonstrated  its  ability  to  perform 
dynamic  plasma  interactions  calculations  for  complex  systems. 

Achievements  in  the  area  of  code  development,  enhancement,  and  validation  include; 

(1)  For  NASCAP/LEO; 

(a)  Demonstration  of  the  ability  to  calculate  equilibrium  sheaths  and  current 
distributions  for  systems  with  high  geometrical  complexity; 

(b)  Addition  of  a  "frozen  ion"  approximation  to  aid  in  analysis  of  transient  sheath 
phenomena. 

(2)  For  Gilbert; 

(a)  Development  of  an  internal  boundary  condition  to  mimic  the  electrostatic  effects 
of  a  screen: 

(b)  Demonstration  of  the  ability  to  predict  spatial  and  temporal  current  distributions 
with  sufficient  accuracy  to  make  meaningful  physical  predictions. 

(3)  For  DynaPAC: 
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(a)  Development  of  a  computer  code  capable  of  performing  dynamic  plasma 
calculations  for  systems  with  high  geometrical  complexity; 

(b)  Development  of  a  programmer-fnendly  DataBase  Manager  for  the  allocation, 
storage,  and  retrieval  of  large  arrays  of  data; 

(c)  Development  of  high-order  finite  element  techniques  to  guarantee  accurate 
potential  solutions  with  strictly  continuous  electric  fields; 

(d)  Development  of  pre-  and  post  processors  for  ease  of  problem  setup,  input 
generation,  and  data  visualization. 

Applications  of  these  codes  to  the  SPEAR  program  include; 

(1)  Prediction  (using  Gilbert)  of  the  rapid  (or  slow)  breakdown  of  a  negatively  biased  high 
voltage  fixture  by  the  presence  (or  absence)  or  calculable  ion  current  to  the  ’triple  point”; 

(2)  Calculation  (using  Gilbert)  of  dynamic  ion  currents  for  the  small  chamber  test  of  SPEAR 
II,  showing  that  the  intermediate-length  pulse  was  most  likely  to  break  down; 

(3)  Calculation  (using  NASCAP/LEO)  of  the  equilibrium  current  distribution  and  estimate  of 
the  equilibrium  floating  potential  for  SPEAR  II; 

(4)  Development  (using  Gilbert)  of  a  theory  for  dynamic  sheath  development  and  transient 
current  for  a  negative  probe,  showing  that  transient  currents  can  exceed  equilibrium 
currents  by  an  order  of  magnitude; 

(5)  Pre-test  prediction  (using  Gilbert  and  NASCAP/LEO)  of  high  parasitic  currents  to  the 
SPEAR  11  high  voltage  probe  in  the  presence  of  plasma; 

(6)  Calculation  (using  DynaPAC)  of  the  time-dependent  plasma  currents  to  the  SPEAR  II  high 
voltage  components. 
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SURFACE  FLASHOVER  THEORY 
WITH  APPLICATIONS  TO  HIGH  VOLTAGE  CONDITIONING 


introduction 

The  SPEAR  program  revealed  the  need  to  condition  high  voltage  components  tor  the  space 
environment.  Both  the  high  voltage  transformer  and  the  Wystrode  bushing  required  a  conditioning 
sequence  before  high  voltage  standoff  could  be  achieved.  Conditioning  is  an  art  and 
unpredictable  and  techniques  that  work  for  one  design  or  environment  can  fail  in  a  different 
situation.  In  practice  it  is  not  possible  to  guarantee  that  a  conditioning  technique  will  work  or  to 
know  how  long  the  conditioning  will  last.  To  gain  insight  into  the  conditioning  process,  a  surface 
flashover  theory  was  developed  that  describes,  quantitatively,  the  initiation  of  fiashover.  The 
theory  reveals  the  role  surface  desorption  plays  in  the  arc  initiation  process  and  is  an  important 
step  in  understanding,  quantitatively  the  process  of  conditioning  model  to  aid  in  the  design  of  high 
voltage  components  and  conditioning  techniques. 

Application 

High  voltage  insulation  in  space  emerges  as  a  critical  issue  because  of  the  advantages  inherent 
in  operating  high  voltage  SDI  systems  in  configurations  where  bushings,  components,  and 
distribution  networks  are  exposed  to  the  ambient  environment.  Recognizing  this,  the  SPEAR 
program  initiated  a  number  of  space  simulation  experiments,  conducted  in  laboratory  plasma 
chambers,  that  demonstrated  that: 

(1)  The  plasma  induced  arcing  could  be  diminished  through  changes  in  the  geometry  of  high 
voltage  components.  Computer  simulations  substantuated  by  laboratory  experiments 
showed  that  the  use  of  electrostatic  screens  and  electrode  shaping  redirected  ions  away 
from  the  triple  points  and  increased  high  voltage  standoff. 

(2)  High  voltage  standoff  capability  of  insulators  improved  with  conditioning.  Conditioning 
consists  of  repeatedly  applying  high  voltages  to  a  component  allowing  it  to  arc  until  it  is 
capable  of  withstanding  the  high  voltage.  The  conditioning  process  is  understood 
qualitatively  as  resulting  from  changes  in  the  surface  properties  such  as  removal  of 
surface  contaminants  and  gas  desorption.  The  theory  developed  for  the  SPEAR  program 
relates  high  voltage  standoff  to  surface  desorption  coefficients  and  conditioning  leakage 
urrents  to  quantity  of  desorbed  gas.  As  it  stands  now  the  theory  can  predict  pressure 
increases  during  the  condition  process  which  can  be  used  as  an  independent 
measurement,  (a)  to  verify  that  conditioning  is  occurring,  and  (b)  to  estimate  the 
conditioning  time  required. 
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Innovative  Concepts 


The  SPEAR  model  of  surface  tiashover  initiation  quantitatively  predicts  breakdown  thresho'cr;.  a.  a 
gas  release  given  surface  electron  impact  deso’ption  properties.  The  theory  combines  th,.-  en  :cis 
of  surface  charging,  surface  secondary  electron  emission,  neutral  desorption,  and  lonizats  s  .•  trim 
the  desorbed  gas.  Because  the  theory  relates  the  quantity  of  desorbed  neutrals  to  the  fas'  over 
current  and  applied  voltage,  it  can  be  tested  in  laboratory  experiments.  With  additiona'  ca'a  cn 
surface  desorption  over  periods  of  electron  Impact,  the  theory  can  be  used  to  predict  tne 
conditioning  needed  to  reduce  the  desorption  rate  below  flashover  thresholds.  This  aisc  cuuia 
be  tested  and  if  confirmed  used  in  the  design  of  high  voltage  components  and  conditioning 
strategies. 

A  Theory  of  Surface  Flashover 

In  SPEAR  chamber  tests  the  neutral  pressure  was  observed  to  increase  during  conditioning  of 
the  high  voltage  components.  This  revealed  the  role  of  surface  neutral  desorption  in  tna  arc 
initiation  process.  Impact  by  impinging  electrons  has  been  proposed  as  the  mechanism  foi  this 
desorption.  We  utilize  the  hypothesis  of  desorption  of  gas  and  its  subsequent  ionization  oy 
electron  bombardment  to  construct  a  model  of  surface  flashover.  The  idealized  configuratic."  in 
Figure  10-18  shows  two  electrodes  separated  by  an  insulator  length  t  and  charged  to  a  potential 
difference  V.  The  inter-eiectrode  space  may  be  filled  by  a  plasma  of  density  Hp. 


See  Figure  10-18  on  page  99 


Figure  10-18.  Idealized  configuration  for  surface  flashover  problem  shows 
insulator  of  length  { with  so.iace  charge  density 

One  imagines  that  the  surface  of  the  insulator  becomes  positively  charged  so  that  any  eiec’xns 
emitted  from  the  cathode  near  the  triple  point  strike  the  insulator,  producing  secondary  eiecti  o  is. 
which  also  remain  close  to  that  surface.  Electrons  bombarding  the  surface  release  adsorbed 
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Figure  10-18.  Idealized  configuration  for  surface  flashover  problem  shows 
insulator  of  length  k  with  surface  charge  density  Og. 
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gases  that  become  ionized  by  collisions  with  sufficiently  energetic  electrons,  leading  to  an 
increase  in  the  flux  of  electrons  bombarding  the  insulator. 


The  model  equations  for  this  system  are  as  follows; 


1  ^ 

e  dt 


Jout  Jill' 


(10-9) 


where  I  is  the  charge  per  unit  area  on  the  insulator  and  and  are  the  ingoing  and  outgoing 
fluxes  of  electrons.  The  fluxes  j,„  and  are  related  by  the  secondary  emission  yield  Y(e)  at  the 
energy  of  the  incident  electrons  is  e. 

U,  ■  Y(.)  j„.  00-10) 

It  is  useful  to  define  the  surface  charge  density,  Og ,  to  be  the  charge  in  a  layer  above  the  surface 

Oe  ""  (10-11) 

where  n,  is  the  electron  density  in  the  space  above  the  layer,  whose  thickness  is  6.  The  charge 
density,  Og,  is  determined  by  continuity  as  shown  in  equation  10-12. 

^  •  f  ■  -i.  •  U  -  s.  (10-12) 

The  three  terms  on  the  right  hand  side  of  equation  10-12  are  the  sources  and  sinks  of  charge  due 
to  charge  leaving  the  layer,  entering  the  layer  from  the  surface,  and  creation  of  charge  in  the  layer 
due  to  ionization.  In  equation  10-12,  K  is  the  flux  y  (cm  ’  sec  ’)  corresponding  to  Og.  The 
ionization  rate  per  unit  area  of  the  current  layer  is  denoted  by  S.  Recombination  and  attachment 
are  neglected. 


For  simplicity  we  have  assumed  that  all  electrons  in  the  layer  are  born  on  the  surface  of  the 
insulator  {y=0)  with  a  single  energy  =  -imv^o  and  a  velocity  v^  in  the  y  direction. 

The  preceding  equation  must  be  complemented  by  equations  that  determine  the  density  N  of  gas 
in  the  layers; 


NV„  •  lo  -  Yl«,  (10-13; 

where  y  is  the  number  of  gas  molecules  desorbed  per  incident  electron  and  Vq  is  the  mean 
velocity  of  desorbed  molecules  normal  to  the  insulator  surface.  It  is  assumed  that  gas  molecules 
move  through  5  in  times  6A/(,  that  are  negligible  in  comparison  with  the  time  scales  of  the 
phenomena  of  interest. 
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We  now  consider  steady  state  solutions  of  the  foregoing  equations  and  seek  conditions  indicative 
of  breakdown.  From  equations  10-9  and  10-10,  we  require 

e  =  e..  (10-14) 


where  e,  is  the  first  crossover  energy  defined  by 

V(e,)  -  1. 


{10-15} 


—  (e,)>0.  (10-16) 

de 


We  can  relate  the  parallel  and  perpendicular  electric  fields  as  follows.  The  secondary  electrons 
leave  the  surface  with  energy  Co  and  return  to  the  surface  with  energy  e,  a  time,  ^  ^ 

eE±  ■ 

The  distance  traveled  under  the  influence  of  the  parallel  electric  filed  is  d  =  EjT  .  This  along  with 
energy  conservation,  e,  -  Cq  =  Ejd,  can  be  used  to  produce 


Ex 


.f  2^0  - 
1^=1  - 


(10-17) 


In  obtaining  this  relation,  one  neglects  electron  energy  loss  by  ionization  as  well  as  by  inelastic 
processes.  The  effect  of  this  neglect  should  be  examined. 

The  ionization  rate,  S,  due  to  collisions  by  the  surface  current  of  electrons,  K,  with  the  desorbed 
gas,  N,  is: 


S  =  NKo  =  -^Ko.  (10-18) 

Vo 

The  surface  current,  K,  can  be  expressed  in  terms  of  the  incident  electron  flux,  J,„  by  noting  that 
K  is  just  the  current  generated  the  secondary  electrons  hopping, 

eK  =  (10-19) 
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2e  E. 


K*  =  aK^. 


llC-20) 


solving  for  K  and  substituting  in  equation  10-19,  we  have 

Substituting  this  into  equation  10-17,  and  gathering  the  constants  into  the  coefficient,  c-,,  gi  es 


ax 


aK' 


(li 


1) 


integrating  equation  10-21  gives 


K  = 


Kc 

1  -  Ko  a  x' 


(10  22) 


where  eKg  is  the  current  at  the  cathode  end  of  the  insulator.  Electrons  may  be  libera'c.3  at  the 
cathode  by  field  emission  or  as  a  result  of  ion  bombardment  of  the  cathode.  In  gerurai.  the 
magnitude  of  Kq  depends  on  the  parameters  of  the  discharge  and  the  surrounding  ambient 
plasma  including,  in  particular,  the  electric  field  and  ion  current  distribution  near  the  cathoce  triDie 
point.  We  denote  these  dependencies  by  the  expression  that  is,  K^,  is  taken  to 

be  some  function  of  the  electric  field  and  ion  current  density  j,  at  the  cathode.  The  latler 
includes  ions  generated  by  collisions  with  electrons  in  the  layer  S  that  move  to  the  cathode  under 
the  influence  of  the  fields. 


Although  the  mechanisms  associated  with  cathode  emission  are  not  well  understood,  we  can 
nevertheless  define  a  critical  emission  current 


corresponding  to  breakdown  conditions. 

Application  to  SPEAR  fl  and  other  Syctems 

We  seek  to  determine  if  the  critical  current  given  in  equations  10-23  is  in  accord  with  the  c  nt 
I  =  1,3  X  10^  amp  observed  by  Gray  in  a  simple  flashover  channel  with  a  gap  f  =  1  3  :Q‘ 

cm,  a  flashover  field  of  20  kV/cm  and  a  channel  width  W'-IO^cm.  This  gives 

qK  = -i  =  .1.^.  (10  24) 

W  cm 

From  equations  15-20 

a,  .  idsl  (VC 

Vo  e,  -  Co 
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Taking  {o)=10'’®  cm^  Vq  =  4  x  10^  cm/sec,  c,  -  cq  =  100  eV,  we  find, 


qK 


1.6  X  10-''^ 
a{ 


JAxJp-^^  .  M  amfVcm. 

1.3  X  10"^®Y  Y 


(10-26) 


This  value  is  in  reasonable  accord  with  experimental  results  for  y~1.  A  desorption  probability  of 
a  few  molecules  per  electron  is  consistent  with  results  reported  by  Hackam  and  is  a 
reasonable  value  if  the  surface  is  heated  by  electrons. 


The  theory  relates  the  surface  current.  K,  to  the  rate  of  desorption  of  gas  from  the  insulator 
surface.  As  will  be  shown  below,  the  rates  predicted  are  in  qualitative  agreement  with  the 
pressure  increases  observed  in  SPEAR  II  during  conditioning  sequences  on  the  pad  at  White 
Sands. 

Using  equation  10-13.  we  can  integrate  the  Hux  of  desorbed  neutrals  over  the  surface  contributing 
to  the  flashover  to  get  the  total  current  of  neutrals,  R. 

R  =  //NVdA^  (10-27) 

surface 

with  equations  10-19  and  10-22  this  can  be  expressed  in  terms  of  the  flashover  current,  1^,  as 

R  =  Y-!^  C*  ~  Kog^)^h(1  -  Kptt^)  (10-28) 

e  Kottd 

During  flashover  described  by  equations  10-22  and  10-23  a  large  but  obviously  finite  current 
flows.  Power  supplies  have  current  limits  and  will  control  K^,  to  give  a  large  but  finite  currents. 
With  this  in  mind  we  can  approximate  equation  10-28  by  assuming  Kat  is  near  1  say  .9-. 99 
corresponding  to  an  amplification  of  10-100.  in  this  case  we  have, 

H  »  (.1  to  .3)  ■  Y 

e  a 


This  equation  for  low  value  of  KgOf  where  the  prefactor  in  equation  10-29  is  one,  states  that  the 


rate  of  desorbed  neutrals  equals  R  = 


Current  of  electrons  )  ( 

in  flashover, -5 

\  e  1  \ 


Number  of  hops 
along  flashover  paths, 


•  Y  The 


total  amount  of  neutrals  produced  during  conditioning  flashover  is  the  rate  of  neutral  production 
times  the  flashover  time.  For  the  SPEAR  II  conditions,  the  number  of  electrons  involved  in  a 
typical  flashover  of  length  x  =  lO'^s  'S 
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rhe  ratio 


^  _  1 03  •  1 0  _  ■JO’* 

©  1.6  •  10  ’®C 


.  -- .  is  the  number  of  hops  of  the  electrons  along  the  (lashover  surface.  T  'us  is 
number  of  intervals  with  energy  e,  50©V  make  up  the  standoff  voltage. 


1  - 

d  '  50v 


200c 


For  each  electron  impact  onto  the  surface  yneutrats  are  released  into  the  chamber.  T?ie  p'cc  se 
value  of  7  is  uncertain  but  is  between  1  and  10  for  impact  energies  of  e,  =  50eV.  The  number 
of  gas  molecules  released  during  a  discharge  is  the  n  =  10'*  -  10'*. 


Assuming  during  the  SPEAR  ii  conditioning  sequences  that  this  amount  of  neuvais  was  emitteo 
into  a  volume  of  about  a  cubic  nteter  it  would  produce  a  pressure  of  between  -10'*  and  10  '  lorr. 
This  is  consistent  with  the  obsen/ations  from  White  Sands  that  the  pressure  increased  cjr<rg 
each  conditioning  sequence. 


Equations  10'20,  10-22  and  equoiion  10-29  are  im.portant  steps  in  understanding  conditioning. 
Equations  10-20  and  10-22  reveal  the  drivers  in  the  surface  flashover  process.  In  general  ;t  s 
not  possible  to  control  the  flashover  path  length.  1,  so  conditioning  must  modify  the  surface 
desorption  properties  described  by  o,  namely  y.  Decreasing  the  surface  desorption  coefficient 
will  increase  tfie  flashover  threshold.  Conditioning  does  exactly  this  by  driving  neutrals  from  the 
Surface  of  the  insulator.  Moreover,  equation  10  29  relates  the  amount  of  surface  desorption  to 
tne  conditioning  current.  With  further  deveioprnent  and  experimental  confirmation,  this  could  be 
used  to  diagnose  conditioning  of  high  voltage  componer.ts  b;  giving  a  measure  to  the  ccndifioning 
process.  For  example,  measuring  the  flashover  current  and  the  pressure  of  the  desorbed  ne  .‘ra'S 
could  be  used  with  equation  10-29  to  give  y  which  could  be  used  as  a  measure  of  the  condition 
of  the  surface. 


New  Technical  Achievements  and  Capabilities 

(1)  A  new  surface  flashover  theory  has  been  developed  for  use  oy  the  SPEAR  program  ana 
other  high  voltage  operations  in  space. 

(2)  The  theory  relates  surface  flashover  current  to  electron  impact  induced  surface  deserptnn 
The  pressure  increases  due  to  these  neutrals  and  are  consistent  with  the  pressure 
observations  during  the  SPEAR  II  conditioning  sequences  on  the  pad  at  White  Sanqs 

(3)  The  theory  provides  an  important  first  step  in  understanding  the  process  of  high  voltage 
conditioning.  Further  development  could  result  in  a  predictive  model  of  conditioning 
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QUICK  MODEL  OF  SPACECRAFT  OUTGASSING  PLUMES 


Introduction 

A  new,  quick  running  model  of  the  bulk  neutral  density  from  outgassing  has  been  developed  after 
the  SPEAR  i  flight  for  the  SPEAR  II  program.  The  model  includes  the  effects  of  reflection  from 
surfaces  as  well  as  scattering  by  the  ambient  ram  neutrals.  The  model  is  quick  running  and  can 
be  integrated  into  engineering  design  codes  such  as  EPSAT‘’*'.  The  model  showed  that  the 
outgassing  of  the  G-10  booms  was  insufficient  to  cause  ionization  breakdown  during  the  SPEAR 
!  flight  even  for  the  extreme  case  where  the  G-10  was  assumed  to  be  responsible  for  the  high 
pressures  measured  by  the  neutral  pressure  gauge.  Additionally,  return  flux  of  outgassing  due 
to  collisions  with  the  ambient  environment  was  shown  to  be  negligible.  Preflight  calculations 
performed  for  SPEAR  II  showed  that  sheath  ionizations  would  have  no  effect  on  the  power 
system  performance.  This  was  confirmed  in  the  Plumbrook  chamber  tests  where  no  sheath 
breakdowns  were  observed. 

Application 

Typical  space  systems  must  operate  in  a  wide  range  of  neutral  contaminant  environments.  The 
ambient  neutral  density  changes  with  season,  time  of  day,  and  orbit.  Attitude  control  thrusters 
periodicaily  fire.  Outgassing  rates  from  surfaces  change  with  time,  and  accommodation  rates 
depend  on  incident  fluxes,  which  in  turn  depend  on  all  of  the  above.  The  neutral  environment  can 
cause  deleterious  effects  to  exposed  high  voltage  components.  Paschen  breakdown  can  occur 
if  the  neutral  density  falls  near  the  Paschen  minimum.  Neutrals  in  electron  colleciing  sheaths  can 
have  bulk  breakdown  (sheath  ionization)  to  the  background  plasma  and  other  components  as  was 
observed  in  the  SPEAR  I  Plumbrook  chamber  tests.  When  designing  for  these  changing 
environments,  it  is  necessary  to  have  quick  estimates  for  the  neutral  density  at  select  positions 
(such  as  near  a  high  voltage  component)  throughout  the  time  period.  The  SPEAR  outgassing 
model  allows  the  rapid  calculation  of  the  bulk  neutral  densities  around  complex  spacecraft  such 
as  SPEAR  i  and  SPEAR  li.  The  model  is  complementary  to  the  Auburn  outgassing  model;  the 
S-Cubed  model  uses  the  Auburn  surface  and  complex  object  rates  as  input  and  calculates  the 
external  expansion  including  scattering  off  surfaces  and  the  ambient  ram. 
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Innovative  Concepts 


The  SPEAR  contamination  model  includes  the  return  flux  scattered  from  the  ambient  ram  The 
model  makes  use  of  the  scale  in  variance  of  the  physics  to  enable  a  quick  look-up  sch>  r  e  to  be 
employed  that  replaces  the  volume  integral  used  in  standard  contamination  codes. 

The  surface  to  surfanA  view  factors,  which  determine  the  reflection  contribution  frcrr,  cije  to 
neutral  accommodation,  are  computed  by  breaking  the  surfaces  into  triangles  and  then  tne 
analytic  result  for  the  solid  angle  subtended  by  a  triangle.  This  approach  replaces  o  s  .rtace 
integral  that  is  used  in  many  contamination  models. 

The  Quick  Outgassing  Model 

The  neutral  effluent  due  to  outgassing  and  accommodation  is  generally  of  low  enough  dor.  ray  mat 
coliisional  mean  free  paths  are  large  compared  with  object  dimensions.  For  this  case,  neutrals 
leaving  a  surface  travel  in  straight  lines  and  either  leave  the  region  of  interest  or  hit  another 
surface.  Neutrals  that  hit  another  surface  either  stick  or  are  accommodated.  This  process  is 
shown  below  in  Figure  10-19. 


See  Figure  lO-l'-J  on  page  107 


Figure  10-19.  The  neutral  density  at  a  point  arises  from  a)  direct  surface  outgassin'-  arh  n) 
reflections  (accommodation)  from  surfaces,  and  c)  scattering  off  of  thu  anb.t  nt 
ram  neutrals. 

The  neutral  density  at  x  is  the  sum  of  the  outgassing  density  from  each  surface  plus  contn&uuors 
due  to  accommodation  of  neutrals  incident  from  other  surfaces  and  the  ram  scattenr.g  E  ach 
surface  contributes  a  density  given  by, 


P  =  //dS  (10  32) 

where  C  determines  the  outgassing  rate.  The  above  expression  is  proportional  to  the  vio/.  hactor 
of  the  surface  S  at  r.  It  is  convenient  to  evaluate  r  directly  above  the  surface,  S 
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Figure  10-19.  The  neutral  density  at  a  point  arises  from  a)  direct  surface 

outgassing,  and  b)  reflections  (accommodation)  from  surfaces, 
and  c)  scattering  off  of  the  ambient  ram  neutrals. 
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p(0)  =  2rX. 


IW33) 


Hence  we  can  express  the  density  at  r  from  S  in  terms  of  the  outgassing  density  at  the  surface 
and  the  solid  angle  Q,  subtended  by  S  at  r. 

p  =  p(0)^.  (10-34) 


The  total  density  from  all  surfaces  is  then, 

p  -i; 

surfaces  I 


{10-35) 


In  the  above  expression,  p,(0)  is  not  the  bare  density  due  to  outgassing  by  surface  t  but  must 
include  accommodated  neutrals  from  other  surfaces.  Designating  p°  to  be  the  bare  outgassing 
density  column  vector  for  all  surfaces  and  p(0)  to  be  a  vector  containing  the  total  neutral 
densities  for  ail  surfaces,  including  multiple  accommodations,  we  have  the  following  nathx 
equation 


P(0) 


=  P*’ 


M-p(0), 


(10  36) 


where  M  is  a  matrix  of  surface  to  surface  view  factors  including  an  accommodation  factor  (which 
conserves  flux  at  the  surface  and  depends  on  the  incoming  flux  and  the  surface  temperature). 
Although  this  can  be  directly  solved  for  p(0)  in  terms  of  p°  as, 

p{0)  =  (^0-37) 

it  is  more  convenient  to  iteratively  estimate  p{0)  using, 

p{0)  =  p°  -  M-p°  *  +...  (10-38) 

Each  iteration  corresponds  to  a  reflection.  Typically  one  iteration  (corresponding  to  single 
reflection)  is  sufficient.  The  calculations  presented  below  use  the  single  reflection  approximation. 
In  practice  M,  is  computed  once  for  a  system  configuration  and  then  used  for  any  ambient  ram 
environments  and  outgassing  rates. 

Finally,  the  total  neutral  density  at  a  point  is  given  by  the  sum  of  the  contributions  from  each 
surface  plus  the  ram  scattering  contribution.  This  latter  component  can  be  expressed  as  a  volume 
integral  over  the  scattering  sources  throughout  space  with 

Pretum  ’  Poutgasslng  3  dOOS0,  (10  39) 
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where  ^  is  the  area  of  the  surface,  Po^j^assmu  'S  the  outgassing  density  at  the 

surface,  k  '  is  apg,„i3„„„  and  V,  ,  V,  are  unit  vectors  in  the  direction  of  the  outgassing  surface  and 
the  ram  respectively.  All  distances  are  expressed  in  terms  of  x,,  the  distance  between  the  surface 
and  observation  point.  In  the  above  equation,  if  the  dot  products  are  less  than  zero,  the  dot 
product  is  taken  to  be  zero.  This  occurs  when  the  integration  volume  is  behind  the  outgassing 
surface  and  when  the  scattering  by  the  ram  cannot  get  to  the  observer. 


Sie  Figure  10-20  on  page 
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Figure  10-20.  The  ram  scattering  contribution  to  the  neutral  density  can  be  expressed  as  a 
volume  integral  over  the  scattering  sources  throughout  space. 

As  seen  in  equation  10-39  the  dimensional  variables  have  been  entirely  scaled  out  leaving  only 
the  orientational  (angular)  dependences  a,  p,  and  tp.  Therefore,  the  integral  can  be  done 
numerically  for  several  values  of  these  angles  and  tabulated  for  subsequent  use.  This  reduction 
of  a  three  dimensional  integral  to  a  simple  table  look-up  is  the  key  to  the  speed  of  the  backscatter 
return  flux  calculations  in  the  SPEAR  outgassing  model. 

Equation  1 0-39  also  shows  that  return  flux  from  outgassing  scattering  is  negligible  unless  the  size 
of  the  source  is  comparable  to  the  mean  free  path  of  the  ambient  neutral.  For  SPEAR  above 
200Km  where  the  neutral  density  is  <10’®/m^ ,  this  is  less  than  1%  effective. 

Application  Tc  SPEAR 

The  SPEAR  experiments  were  designed  to  investigate  the  space  environment  interactions  with 
high  voltage  components.  In  SPEAR  I,  high  voltage  plasma  current  collection  in  the  presence  of 
the  earth's  magnetic  field  was  investigated.  At  issue  were  the  extent  that  non-classical 
mechanisms  such  as  anomalous  scattering,  turbulence,  and  ionization  affect  current  collection. 
The  methodology  of  this  investigation  was  to  compare  the  SPEAR  I  experimental  results  for 
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Figure  10-20.  The  ram  scattering  contribution  to  the  neutral  density  can  be 
expressed  as  a  volume  integral  over  the  scattering  sources 
throughout  space. 
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current  collection  with  the  results  of  the  NASCAP/LEO^"^^  and  POLAR'^^'  codes  which  assume 
classical  current  collection  mechanisms  (plasma  particles  move  under  the  influence  of 
non-oscillating  fields).  This,  however,  would  be  conclusive  only  if  the  effects  of  neutral  ionization 
could  be  ignored.  The  neutral  density  calculations  using  the  outgassing  model  discussed  showed 
just  that.  Preflight  calculations  were  performed  for  SPEAR  II.  The  calculations  showed  that 
external  outgassing  would  not  cause  sheath  ionization  breakdown. 

Figure  10-21  shows  the  SPEAR  I  object  and  ranges  for  surface  outgassing  rates  provided  by 
Lloyd  Gordon  of  Auburn  University  for  the  hot  aluminum  body,  G-10  boom  and  G-10.'aluminum 
bushings  and  aluminum  spheres.  The  major  unknown  in  these  calculations  was  the  outgassing 
rate  of  the  G-10  booms.  Figure  10-22  shows  the  neutral  densities  about  SPEAR  I  as  calculated 
using  the  above  theory.  The  highest  densities  in  this  calculation  are  less  than  I0’®m  ^  This  is  far 
lower  than  the  neutral  densities  measured  during  the  SPEAR  I  flight  (>10'^m^).  Because  the 
G-10  outgassing  rate  was  suspect,  we  increased  the  G-10  rate  until  the  neutral  density  in  the 
vicinity  of  the  neutral  gauge  reached  the  observed  levels.  This  required  an  increase  in  the  G-10 
outgassing  rate  by  two  orders  of  magnitude.  This  new  outgassing  rate  was  used  in  the  calculation 
shown  in  Figure  10-23.  This  calculation  was  considered  an  extreme  upper  bound  on  the  neutral 
density  about  SPEAR  I.  As  seen  in  the  figure,  even  for  the  extreme  outgassing  rate  assumed  for 
the  G-10  boom,  in  the  sphere  region  the  neutral  density  is  significantly  less  than  10’^m^  For 
these  densities  ionization  effects  in  electron  sheaths  are  not  important. 

In  summary,  outgassing  on  SPEAR  I  could  not  have  caused  any  anomalous  plasma  effects.  This 
is  true  even  for  the  extreme  assumption  that  the  g-10  boom  was  highly  outgassing  and 
responsible  for  the  high  neutral  pressure  gauge  readings  during  flight. 

Calculations  of  SPEAR  II  external  outgassing  densities  were  performed  assuming  the  same 
outgassing  rate  for  ail  materials.  The  value  chosen  was  0.01  W/m^,  the  maximum  ambient  rate 
for  hot  aluminum.  The  resultant  gas  cloud  was  orders  of  magnitude  less  dense  than  necessary 
for  ionization  enhanced  sheath  currents  of  for  sheath  breakdown.  The  calculated  gas  cloud 
surrounding  SPEAR  II  is  shown  in  Figure  10-24.  The  SPEAR  II  Plumbrook  chamber  tests  showed 
no  bulk  ionization  breakdown  in  the  sheath,  consistent  with  the  low  outgassing  densities  predicted 
by  the  model. 


Ill 


Units  X  10^W/m^ 

Aluminum  sphere  with  trapped  gas  5,000  -  7,000 


Figure  10*21.  Outgassing  rates  for  SPEAR  1  components. 
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Figure  10-22.  Outgassing  densities  about  SPEAR  I,  taking  outgassing  rates  of 
7  X  1 0“^  W/m^  for  the  sp!  leres,  1 0'^  W/m^  for  the  bushings, 

5  -  1 0’^  W/m^  for  the  booms,  and  1 0'^  W/m^  for  the  body.  These 
values  were  provided  by  Auburn  University. 
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Figure  10-23.  Outgassing  densities  about  SPEAR  I  using  elevated  G-10 

outgassing.  Even  for  this  enhanced  outgassing  and  resultant  neutral 
densities,  the  density  in  the  sphere  region  is  below  the  sheath 
ionization  thresholds. 
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See  Figure  10-24  on  page  116 


Figure  10-24.  Outgassing  densities  about  SPEAR  II  using  outgassing  rates  of  .01  W/m^  for  all 
surfaces. 

Summary 

The  new  SPEAR  mode!  enables  the  neutral  density  due  to  outgassing,  including  scattering  from 
the  ambient,  to  be  computed  quickly  making  system  trade  studies  practical.  The  model  has  been 
incorporated  into  the  NASCAP/LEO  and  EPSAT  codes  for  use  in  engineering  design  analysis. 
Calculations  using  the  model  show  that  normal  outgassing  has  little  impact  on  collecting  sheaths 
for  the  power  systems,  and  that  ambinet  neutral-outgassing  scattering  is  negligible  for  few  meter 
size  systems  above  200  km.  For  SPEAR  I  the  neutral  density  in  the  electron  collecting  sheath 
due  to  outgassing  was  shown  to  be  below  breakdown  thresholds.  This  result  is  consistent  with 
the  lack  of  sheath  ionization  breakdown  during  hte  SPEAR  1  flight.  For  SPEAR  II,  the  neutral 
density  in  the  ion  collecting  sheaths  was  shown  to  be  insufficient  to  cause  significant  ionization 
enhancements  to  the  collected  current.  This  is  consistent  with  the  Plumbook  chamber  tests.  In 
Plumbrook  the  observed  current  collection  was  in  agreement  with  DynaPAC  calculations,  implying 
that  only  classical  phenomena,  not  ionization,  controlled  the  plasma  interaction  physics. 
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Figure  10-24.  Outgassing  dens 
.01  W/m^  for  all  s 
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Comparison  Matrix 
NASCAP/LEO  -  Gilbert  -  DynaPAC 


NASCAP/LEO 

Gilben 

DynaPAC 

Dimensionality 

3-D 

2V2-D 

3-D 

Potential  Interpolation 

TriLinear 

BiQuadraiic 

High  Order 

Steady-State  Sheaths 

Yes 

Yes 

Yes  1 

Generalized  Geometry 

Yes 

Yes 

Yes  j 

1 

Time  Dependence 

No 

Yes 

\’es  ! 

Magnetic  Fields 

Yes 

ves 

Yes 

Self-Consistent  Sheaths 

No 

f" - -  ■■■  ' 

Yes 

Yes 

Plasma  Wakes 

Yes 

No 

Yes 

Neutral  Ionization 

No 

Yes 

Planned 

DataBase  Centered 

No 

No 

Yes 

User  Intert'ace 

Good 

Good 

Best 

i 

Development  for  SPEAR  II 

Augmented 

Augmented 

Developed  i 

Validation  by  SPEAR  II 

No 

Yes 

Yes 

Extend  SPEAR  11  Results 

Yes 

Yes 

Yes 

Note  to  Editor: 

This  table  is  provided  for  editorial  convenience.  It  does  not  exist  on  the  disks, 
and  is  not  called  out  in  the  text  It  may  be  conveniently  placed  in  the  '’Summar>’ 
of  Technical  Developments  and  Capabilities"  section  at  the  end  of  the  code  por¬ 
tion  of  the  chapter. 
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AppifM  High  VoJtage 
anci 

:  etectronics 


8.  SIMPLE  CALCULATIONS  FOR  THE  BREAKDOWN  OF  A  NEGATIVE 

SHEATH 


This  work  appeared  in  the  Quarterly  Report  for  1  January  through  31  March 
1991. 
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8.  IONIZATION  CALCULATIONS 


Ionization  breakdown  for  negative  potentials  fakes  place  according  to  the 
following  mechanism; 

1 .  An  electron  emitted  from  the  probe  surface  creates  r  ions  as  it  is 
accelerated  into  space; 

2,  The  collection  of  those  ions  results  in  emission  of  electrons. 

We  have  breakdown  for  7r>  1.  More  precisely,  the  current  amplification  factor 
for  a  “seed  ion’’  is  1  +  yo  (f  +  r  )/(l  -  yO.  where  yo 's  the  secondary  yield  for  the 
seed  ion.  Note  that: 

1 .  r  is  a  function  of  the  potential  distribution  V{r),  the  neutral  density 
distribution  Nof(r),  and  the  ionization  cross-section  o(E). 

2.  Y  'S  a  function  of  the  distribution  of  potentials  at  which  ions  are  created. 

One-dimensional  calculations  have  been  done  under  the  following 
assumptions; 

1.  The  probe  is  represented  by  a  sphere  of  radius  R  =  0.64  meters  (hopefully, 
an  adequate  representation  of  the  SPEAR  I  rocket  body). 

2.  The  ionization  cross-section  is  that  for  electron  impact  ionization  of  N2- 

3.  The  potential  is  a  confined  coulomb  potential; 

V(r)  =  Vo  (RfiLB/r  ■  R)/(Rlb  *  R) 

where  Rlb  is  the  Langmuir-Blodgett  radius  for  a  plasma  of  density  10^'' 
and  temperature  0.1  eV. 

4.  The  neutral  density  function  is  f(r)  =  1  or  f(r)  =  (R/r)2. 
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5.  The  secondary  yield  is 


Yo(v)  =  Max[0-05.  2.4  ^  10-5(v  -  5.5  >  10^)] 


where  v  is  the  ion  incident  velocity. 


6.  tons  and  electrons  are  created  with  zero  energy  and  with  no  energy  loss  to 
the  primary. 


The  results  turn  out  to  be  only  a  little  bit  nonlinear.  Here  they  are; 


;  Vo 

Rlb 

fir) 

Y 

r 

Breakdown  Nq  i 

'  -1000 

3  3  1 

1 

0  84+7  5x10-21No 

3.6x10-20No+1  1x10-39No2 

17x10^8 

;  -10. coo 

8  1 

1 

5.68+4  5x10-2ONq 

1  9x10'20No+1  1x10'39No2 

7x10’8 

I  -50,000 

15.6: 

1 

15  3+1  OxIO-’^No 

1.1x10-20Nc+8.9x10-'^°No2 

4x10^8  i 

i  -1000 

3.3 

(R/r)2 

0  435+7.5x10-21No 

9.3x10-21No+6.7x10-^’No2 

1  0x1020  i 

'  -10,000 

8.1  j 

{R'r)2 

2.72+4.5x10-2ONo 

3. 1x10-2’ No+8  8x10-‘^2nq2 

0  8x1020  i 

i  -50,000 

15.61 

(R/r)2 

7.44  +  1. OxIO'^^No 

1.0x10-21No+1.2x10-''2nq2 

1.1x1020  1 

A  few  additional  comments; 

1 .  For  the  {R/r)^  density,  all  the  ion  production  takes  place  within  about  two 
meters  of  the  probe  (except  for  the  -1000  volt  case). 

2.  For  constant  density,  as  breakdown  is  approached,  we  get  increasing 
ionization  at  the  outer  boundary. 

3.  When  we  exceed  about  half  the  breakdown  density,  we  get  order-of- 
magnitude  current  enhancement, 

4.  Breakdown  develops  on  an  ion  transit  timescale,  which  would  be  tens  of 
microseconds.  Accumulated  ions  will  sharpen  the  potential  during  this 
process,  reducing  r  from  what  we  have  calculated,  and  leading  to  a  higher 
breakdown  density. 

5.  In  positive  potential  breakdown,  the  accumulated  ions  force  trapping  of 
additional  electrons  to  balance  their  charge.  For  negative  potential 
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breakdown  at  threshold,  the  electron  density  is  negligible  and  the  retained 
ions  tend  to  suppress  breakdown  by  causing  electrons  to  be  accelerated 
beyond  the  energy  for  maximum  cross-section  more  quickly. 

6.  Calculations  can  proceed  on  an  ion  timescale  by  generating  the  entire 
ionization  track  of  a  secondary  electron  each  time  an  ion  strikes  the  probe 
This  approach  is  valid  as  long  as 

(1/Q)(dQ/dt) »  1/Te 

where  Q  is  the  ion  content  of  the  grid,  and  tg  is  the  time  for  an  ionization- 
produced  electron  (not  a  secondary  electron)  to  leave  the  grid.  That  is,  the 
density  must  be  near  the  breakdown  threshold  rather  than  well  above  it 
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9.  SPEAR-3  SIIEA7  fi  S TRUCrifKE  A\D  PLASMA  C  URREN  I  S 


T!iis  \\.t>rk  iii  tlu*  Qinirtcriv  RrpDrt  lor  I  April  llirouj^E  3il  juno  PM  1 


9.  SPEAR  ill  CALCULATIONS  FOR  PDR 


DynaPAC  was  used  to  calculate  sheath  structure  and  plasma  currents  for  the 
SPEAR  III  preliminary  design  review.  Logan,  Utah.  12-14  June,  1991.  A 
PATRAN  model  was  constructed  having  the  correct  rocket  dimensions.  Half  the 
gradings  were  removed  from  the  boom,  as  had  been  suggested  by  the 
mechanical  engineers.  As  the  boom  angle  was  not  decided  at  the  kickoff 
meeting,  it  was  set  to  approximately  45  degrees.  The  model  is  shown  in 
figure  9.1. 

Steady-state  calculations  were  done  for  sphere  biases  of  5  kilovolts  and  20 
kilovolts  and  several  rocket  potentials  for  each  case.  The  plasma  was  taken  as 
0+  with  density  1  x  IC^  m-3  and  temperature  0.1  eV  and  the  magnetic  field  was 
0.4  gauss.  Figures  9.2  show  potentials  calculated  for  20  kilovolt  sphere  bias, 
and  figures  9.3  show  potentials  calculated  for  5  kilovolt  sphere  bias. 

Tables  9.1  and  9.2  show  the  calculated  electron  and  ion  currents.  Ion  currents 
were  calculated  using  the  usual  sheath  method,  and  no  problems  were 
encountered  except  for  the  two  cases  when  the  sheath  filled  the  grid.  Figure  9.4 
shows  ion  trajectories  for  the  case  of  20  kilovolt  bias  and  negative  4  kilovolt 
rocket  potential.  We  also  added  the  ObjPotl  the  capability  of  plotting  cell 
currents  and  surface  current  densities.  These  are  shown  (for  the  same  case)  in 
figures  9  5.  Peak  ion  current  of  about  0.7  milliamperes  per  square  meter  occurs 
near  the  middle  of  the  rocket  and  opposite  the  sphere,  .  he  high-voltage  section 
(top  section)  of  the  rocket  is  shielded  from  ion  current  by  the  positive  sheath, 
particularly  on  the  side  toward  the  sphere. 
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Table  9.1  Current  (milliamperes)  results  for  SPEAR  IN  with 
magnetic  field  normal  to  rocket-boom  plane. 


Sphere  Body  |  Electron  Current 

I  Fieldline  Method 

Potential  Potential  Sphere  Boom  Tra 


19.5  kV 


18.0  kV 


16.0  kV 


14.0  kV 


12.0  kV  -8.0  kV  I 


5000  V 


Ion  Secondary  Total 
Current  Electrons 


4000  V  -1000  V 


2500  V I  •2500  V I 
2000  V I  -3000  V 


Table  9.2  Current  (milliamperes)  results  for  SPEAR  Iff 
with  magnetic  I'ield  along  boom. 


Ion  Secondary  i  Total 
Current  Electrons  | 


16.0  kV 


14.0  kV 


12.0  kV 


5000  V 


4500  V 


4000  V 


3500 

_ 3000 

2500 

2000 


-6.0  kV 


-8.0  kV 


-500  V 


1000  V 


-1500  V 


-2000  V 


Electron  currents  to  the  sphere  and  boom  proved  difficult  to  calculate.  The 
sheath  method  gave  implausible  results  and  trajectories,  particularly  for  cases 
where  the  potential  was  not  extremely  well  converged.  This  is  due  in  part  to 
unphysical  initial  conditions  given  electrons  started  by  the  sheath  method. 
Eventually,  we  elected  to  follow  electrons  along  field  lines  from  the 
computational  boundary.  This  had  the  disadvantage  in  that  a  lot  of  effort  went 
into  tracking  electrons  through  electric-field-free  regions  but  had  the  nice 
properties  that  trapped,  as  well  as  collected,  currents  that  could  be  calculated, 
and  they  did  not  depend  (within  reason)  on  the  initial  particle  energy. 
(Electrons  were  considered  quasi-trapp  d  if  they  survived  after  three  grid- 
index  loops.)  Figure  9.6  shows  trajectc  --s  for  electrons  emitted  along  the  line 
Z=1,  X=8.5  and  shows  examples  of  unperturbed,  quasi-trapped,  and  collected 
electrons. 

The  electron  c..'  is,  as  expected,  were  very  different  between  the  cases  that 
the  magne»'o  field  was  (1 )  normal  to  the  boom-rocket  plane  or  (2)  along  the 
boom.  I:.  the  first  case,  collected  currents  were  somewhat  in  excess  of  the 
Parkc^r-Murphy  limit,  with  comparable  currents  collected  by  the  boom  and  the 
sphere,  and  the  quasi-trapped  current  was  relatively  small.  In  the  second 
case,  the  sphere  collected  well  under  the  Parker-Murphy  limit,  the  boom 
collected  negligible  current,  and  the  quasi-trapped  current  was  very  high. 
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Fiaure  9.2  SPEAR  111  sheath  contours  for  20  kilovolt  bias  and  rocket  body  at 
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(a) 

Figure  9.3  SPEAR  III  sheath  contours  for  5  kilovolt  bias  and  rocket  body  at  (a) 
zero;  (b)  -500  volts;  (c)  -1000  volts;  (d)  -1500  volts. 
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Figure  9.5  ObjPotl  sutlace  ion  current  plots  tor  case  of  figure  9.4  and  9.2(c); 
(a)  surface  cell  currents;  (b)  surface  current  densities. 


Spear3  +16000  -4000 


COLOR  LEGEND 
8  E-04  __ 


7  5E-04 , 
7E04 
6  5E-04 
6  6-04 
5  5E  04 
5.E04 
4.SE04 
4  E  04 
3.5E-04 
3  E-04 
2.5E-04 
2.E-04 
1.5E-04 
l.E-04 
5.E05 
E-»00 
-5.E05 


(b) 


140 


Y-AXIS 


15. 


Siic9  Z  -  8.CQQQ 

Mesh  Size  -  5.0000E-01  meters 


Date-  05-31-91 
Time:  1 1 :43:03 

cot OR  LEGEND 


Figure  9.6  Electron  trajectories  for  potential  of  figure  9.2(c)  and  magnetic  field 
of  0.4  gauss  along  Z  axis.  Electrons  are  emitted  along  field  lines 
from  computational  boundary  at  Z=1 ,  X=8.5. 


10.  SPEAR-3  CALCULATIONS  FOR  PARTICLE  DETECTOR  PLACEMENT 


This  work  appeared  in  the  Quarterly  Report  for  1  April  through  30  June  1991. 
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10.  SPEAR  ill  CALCULATIONS  FOR  PARTICLE 
DETECTOR  PLACEMENT 


Following  the  SPEAR  III  PDR,  we  were  requested  to  do  calculations  to  assist  in 
determining  the  placement,  view  angle,  and  angular  aperture  of  the  ion 
detectors.  As  the  strawman  configuration  had  changed  to  a  right-angle  boom,  a 
new  model  was  generated  and  gridding  was  improved  to  be  suitable  for  any 
bias  in  the  range  of  interest. 

We  performed  calculations  for  the  angular  distribution  of  incident  ions  for  three 
representative  sheath  configurations.  We  recommended  that  if  the  ion 
detector(s)  are  located  on  Science- 1 ,  that  they  be  located  on  the  side  away  from 
the  sphere  and  be  aimed  upward  to  a  degree  depending  on  their  precise 
location. 

The  three  potential  conditions  were: 

1 .  16  kV  bias,  body  at  -4  kV. 

2.  16  kV  bias,  body  at -500  volts. 

3.  5  kV  bias,  body  at  -2500  volts. 

The  plasma  density  and  temperature  were  2  x  10^^  m*3  and  0.1  eV, 
respectively.  A  packet  of  plots  was  generated  for  each  condition.  Figures  10.1 
refer  to  condition  1 ,  figures  10.2  to  condition  2,  and  figures  10.3  to  condition  3. 

The  calculations  were  done  in  three  ways: 

1 .  The  rocket  was  taken  to  be  a  cylinder  with  the  boom  and  sphere  ignored. 
The  Gilbert  code  was  used  to  make  2-dimensional  (R-Z)  potential  and 
trajectory  calculations.  Results  are  included  as  line  plots  (angle  vs.  inches) 
in  the  figures  for  conditions  (1)  and  (2). 
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2.  Particles  were  tracked  inward  from  the  sheath  to  the  rocket  using 
DynaPAC.  These  calculations  were  done  in  the  plane  of  the  potential  plots. 
Each  packet  contains  a  particle  trajectory  plot  and  an  angle  vs.  inches  plot 
with  circles  and  (barely  visible)  crosses.  The  circles  are  results  for  the  side 
of  the  rocket  opposite  the  sphere  and  the  crosses  for  the  side  toward  the 
sphere.  This  method  is  particularly  useful  in  determining  which  regions  of 
the  rocket  receive  no  incident  ions. 

3.  Particles  were  tracked  outward  from  proposed  detector  locations  12  inches 
from  either  end  of  Science-I.  The  results  are  summarized  in  table  10.1. 
Each  set  of  figures  contains  two  trajectory  plots  for  each  of  the  two 
locations.  Trajectories  that  escape  the  sheath  represent  possible  incident 
ion  trajectories:  trajectories  that  reimpact  the  rocket  do  not.  The  trajectory 
plot  in  the  plane  of  the  potentials  shows  incident  ions  arriving  at  a  narrow 
range  of  angles.  The  centra!  trajectory  was  rotated  about  a  line  through  the 
detector  parallel  to  the  rocket  axis,  showing  that  the  range  of  azimuthal 
angles  is  relatively  large. 


Table  10.1 

DynaPAC  reverse  tracking  results  for  ion  detector  angle.  Detector  locations  are 
on  Science- 1 ,  on  side  of  rocket  opposite  right-angle  boom,  with  location  in 
inches  as  indicated  in  distributed  drawings.  Locations  are  12  inches  above 
bottom  of  Science-1  and  12  inches  below  top  of  Science-1.  Angle  8  is  angle 
from  normal  to  rocket  toward  nose.  Angle  $  is  rotation  from  normal  to  rocket 
about  a  line  on  surface  parallel  to  rocket  axis.  The  ±  angles  are  not  errors  but 
are  our  computational  estimates  of  the  angular  spread  of  the  beam.  They  are  in 
reasonable  agreement  with  our  analytic  estimates. 


Sphere 

Body 

Detector 

e 

0 

Bias 

Potential 

Location 

(deg.) 
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The  three  computational  methods  are  in  good  agreement.  Detectors  on 
Science- 1  should  be  located  on  the  side  of  the  rocket  away  from  the  sphere. 

The  incident  angle  increases  with  distance  from  the  rocket  center  ^'s^gniy  the 
bottom  of  Science  1).  The  rate  of  increase  is  about  0.1  degrees  per  inch  with  the 
rocket  at  -500  volts  and  about  0.3  degrees  per  inch  with  the  rocket  at  -4000 
volts.  An  additional  Gilbert  calculation  (figure  10.4)  indicates  that  the  angle  does 
not  increase  too  much  further  for  a  very  negative  (-12  kV)  rocket.  A  detector  on 
the  sphere  side  will  have  data  dropout  when  the  sphere  is  at  high  potential,  and 
a  90  degree  location  has  a  high  risk  of  data  dropout.  A  15=  aperture  detector 
near  the  bottom  of  Science-1  should  cover  the  range  -5=  to  -t-lO”  upwards.  Near 
the  top  it  should  cover  +5=  to  -t-20=.  Azimuthally,  the  detector  should  be  aimed 
straight  outward.  We  would  expect  such  a  detector  to  receive  essentially  the 
entire  incident  ion  flux  to  its  active  area. 

It  is  interesting  that  the  azimuthal  angle  range  for  the  16  kV  bias,  -500  volt  rocket 
case  appears  anomalously  large.  A  contour  plot  (included  in  figure  10.2)  normal 
to  the  rocket  axis  at  the  level  where  incident  ions  enter  the  sheath  shows  that 
the  influence  of  the  positive  sphere  makes  the  ion  collecting  sheath  more  nearly 
concentric  with  the  detector  location  than  with  the  rocket  axis,  explaining  the 
large  range  of  azimuthal  angles. 
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Figure  10.1(a)  SPEAR  III  potentials  for  right-angle  boom,  20  kilovolt  bias,  -4 
kilovolt  rocket. 
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Figure  10.1(e)  DynaPAC  reverse  trajectory  calculation  for  ion  detector  near 
top  of  Science-1 . 
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Figure  10.1(f)  DynaPAC  reverse  trajectory  calculation  for  azimuthal  angle 
range  corresponding  to  figure  10.1(e). 
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Figure  10.1(h)  DynaPAC  reverse  trajectory  calculation  for  azimuthal  angle 
range  corresponding  to  figure  10.1(g). 
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Figure  10.2(a)  SPEAR  Hi  potentials  for  right-angle  boom,  20  kilovolt  bias,  -500 
volt  rocket. 
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Figure  10.2(f)  DynaPAC  reverse  trajectory  calculation  for  azimuthal  angle 
range  corresponding  to  figure  10.2(e). 
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Figure  10.2(h)  DynaPAC  reverse  trajectory  calculation  for  azimuthal  angle 
range  corresponding  to  figure  10.2(g). 
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Figure  10.2(i)  Potential  contours  in  a  cut  normal  to  the  rocket  axis,  showing 
non-concentric  sheath  which  leads  to  a  wide  range  of 
azimuthal  angles. 
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Figure  10.3(d)  DynaPAC  reverse  trajectory  calcuiation  for  ion  detector  near 
top  of  Science-1 . 
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Figure  10.3(e)  DynaPAC  reverse  trajectory  calculation  for  azimuthal  angle 
range  corresponding  to  figure  10.3(d). 
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Figure  10.3(t)  DynaPAC  reverse  trajectory  calculation  for  ion  detector  near 
bottom  of  Science-1. 
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Figure  10.3(g)  DynaPAC  reverse  trajectory  calculation  for  azimuthal  angle 
range  corresponding  to  figure  10.2(f). 


Y  AXI 


11.  THREE-DIMENSIONAL  MAGNETIC  SHEATH  CALCULATIONS 


We  used  the  DynaPAC  computer  code  to  perform  3-dimensional  calculations  of  a 
0.8  m  radius,  1000  V  sphere  moving  through  a  plasma  with  a  magnetic  field.  The 
orbital  motion  was  taken  as  7500  m/sec  in  the  x  direction,  and  the  magnetic  field 
was  0.4  gauss  in  the  y  direction.  The  plasma  density  was 
10'’2  m-3,  and  the  plasma  temperature  was  0.1  eV.  The  ion  species  was  taken  to 
be  0+. 

Figure  11.1  shows  the  wake  of  the  uncharged  object.  The  calculation  uses  a 
shadowing  algorithm  (originally  developed  for  the  POLAR  code)  to  obtain  the 
“neutral  approximation”  results,  then  applies  an  electric  field  correction. 
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Figure  11.1  Plasma  wake  densities  for  an  uncharged,  0.8  m  radius  sphere. 
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Figure  1 1 .2  shows  the  electrostatic  potential  calculated  using  the  object-wake  ion 
densities.  The  calculation  takes  account  of  reduced  screening  due  to  low  ion 
density  in  the  wake  as  well  as  reduced  wake-side  electron  currents.  Thus,  the 
high  positive  potential  extends  into  the  wake. 
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Figure  11.2  Electrostatic  potentials  calculated  with  the  sphere  at  1000  volts  and 
the  ion  densities  of  Figure  11.1. 

It  is  apparent  from  Figure  1 1 .2  that  the  initial  wake  calculation  is  wrong.  The 
wake  of  a  positive  object  is  formed  not  by  the  object  surface,  but  by  the  much 
larger  ion  stagnation  surface.  The  ion  stagnation  surface  coincides  with  the  5  V 
contour  level  in  the  ram  direction  and  drops  to  lower  potential  values  on  the  sides 
of  the  sheath  since  the  electric  field  is  oblique  to  the  ram  velocity.  Figure  1 1 .3 
shows  wake  ion  densities  calculated  as  in  Figure  11.1  but  now  using  the  ion 
stagnation  surface  as  the  shadowing  surface. 
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Figure  11.3  ion  densities  calculated  using  the  ion  stagnation  surface  as  a 
shadowing  surface. 

Figure  1 1 .4  shows  potentials  calculated  using  the  wake  ion  densities  of 
Figure  1 1 .3.  Also,  the  calculation  uses  a  new  space-charge  formulation  that 
allows  the  wake  region  to  be  electron  rich.  We  see  the  positive  potentials 
extending  further  into  the  wake,  and  negative  potentials  down  to  about  a  half  volt 
negative. 
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Figure  11.4  Electrostatic  potentials  calculated  using  ion  densities  similar  to 

Figure  32,  end  a  space  charge  formulation  allowing  the  wake  to  be 
electron  nch 

Finally,  we  add  a  magnetic  field  along  the  Y  axis.  To  determine  self-consistent 
potentials  in  a  magnetic  field  requires  generating  electrons  at  a  "sheath  surface" 
and  tracking  them  to  determine  seff-consisfent  electron  densities  within  the 
sheath.  Figures  1 1  .Slabo  show  me  resultant  electron  densities,  and 
Figures  1 1 .6(abc)  show  the  potentials  fpe  three  electron  density  plots  show  low 
electron  densities  near  the  -,Y  sides  ot  the  sphere  but  relatively  higher  densities 
in  the  other  directions.  This  indicates  that  a  belt  of  electrons  is  now  orbiting  the 
sphere  in  the  X*Z  plane.  T".  s  belt  of  electrons  has  little  effect  on  the  ram  side 
(where  its  charge  is  n'^utrai.zed  by  t.he  ions),  but  on  ttie  wake  side  (where  there 
were  previously  no  ions  or  olectrcns)  it  squeezes  the  sheath  boundary  close  to 
the  object,  as  is  seari  by  comparing  Figufes  1  l.6{abc)with  Figure  1 1.4.  It  also 
reduces  the  Z  sheath  dimension  ’'elative  to  the  Y  sheath  dimension,  although  for 
these  parameters  the  asymmetry  only  b£..omes  apparent  wakeward  of  the 
sphere.  (Rough  caicuiation?,  for  other  parameters  show  far  more  apparent 
asymmetry.) 
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Figure  11.5(a)  Tracked  electron  Densities  in  the  sheath  region,  with  a  magnetic 
field  of  0.4  gauss  in  the  y  direction:  (a)  X-Y  plane. 
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Figure  11.5(b)  Tracked  electron  densities  in  the  sheath  region,  with  a  magnetic 
field  of  0.4  gauss  m  the  y  direction:  (b)  X-Z  plane. 
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Figure  11.5(c) 


Trnck  .t :  ;rc’'cr.  dynsit'E-s  sn  the  sheath  region,  with  a  magnetic 
fieid  ot  L-  4  gauss  m  the  y  direction;  (c)  Y-Z  plane. 
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Figure  11.6(b)  E' 
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sphere  Wake  C  a;.  ui;-'.,-; 
Slice  X=  29  0000 

Units;  METERS 
Min* ‘8.48 1 4E'CJ 

10, 

9.5 
9, 

8.5 
8, 

7.5 
7. 

6.5 
6. 

5.5 
5. 

4.5 
4, 

3.5 
3, 

2.5 
2, 

1.5 
1, 


Date:  03-24  92 
Time:  16:58:24 


COLOR  LEGEND 


-100. 
-0.5, 
-0.4 
-0.3 
-02 
-0  1 
.1 
1 

5. 

20 

i 

200  1 

Jt/;; 

400. 

600 

800 

1000, 


Figure  11.6(c) 


